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Avant-Propos 



When studying any physical problem in Applied Mathematics, three es- 1 
sential stage are involved. 

1. Modelling: An appropriate mathematical model, based on the 
physics or the engineering of the sitution, must be found. Usu- 
alluy these models are given a pariori by the physicists or the 
engineers themselves. However, mathematicians can also play an 
important role in this process especially considering the increas- 
ing emphasis on non - linear models of physical problems. 

2. Mathematical study of the model: A model usally involves a set of 
ordinary' or partial differential equations or an (energy) functional 
to be minimized. One of the first tasks is to find a suitable func- 
tional space in which to study the problem. Then comes the study 
of existence and uniqueness or non -uniqueness of solutions. An 
important feature of linear theories is the existence of unique so- 
lutions depending continuoussly on the data (Hadamard's defini- 
tion of well - posed problems). But with non-linear problems, 
non-uniqueness is a prevealent phenomenon. For instance, bifu- 
racation of solutions is of special interest. 

3. Numerical analysis of the model: By this is meant the descrip- 
tion of, and the mathematical analysis of, approximation schemes, 
which can be run on a computer in a 'reasonable' time to get 'rea- 
sonably accurate' answers. 

In the following set of lectures the first two of the above aspects will 
be studied with reference to the theory of elasticity in three dimensions. 



v 
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2 In the first chapter a non-linear system of partial differential equa- 

tions will be established as a mathematical model of elasticity. The 
non-linearity will appear in the highest order terms and this is an impor- 
tant source of difficulties. An energy functional will be established and 
it will be seen that the equations of equilibrium can be obtained as the 
Euler equations starting from the energy functional. 

Existence results will be studied in the second chapter. The two 
important tools will be the use of the implicit function theorem and the 
theory of J. BALL. 
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Chapter 1 



Description of Three - 
Dimensional Elasticity 

THIS CHAPTER WILL be divided into four sections. In the first sec- 3 
tion some preliminaries on deformations in K 3 will be discussed; the 
second will be devoted to the equations of eqilibrium and the third to 
constitutive equations. These together will give rise to the boundary 
value problem which will serve as the model for three - dimensional 
elasticity. The last section will describe the energy functional and the 
associated Euler equations will be seen to give the equations of equilin- 
rium and the constitutive equations. 

1.1 Geometrical Preliminaries 

Let O c R 3 be a bounded open set. Let 23# = Cl, the closure of f2 in 
R 3 , stand for the reference configuration. (The subsript R will always 
stand for the reference configuration.) Let Xr be a generic point in 23#. 
If [e\,e2, ej] is the standard orthonormal basis for M 3 , 

(1.1-1) OX R = X R , ei 

where OXr stands for the position vector of Xr. (In the above relation 4 
and in all that follows, the summation convention for repeated indices 
will always be adopted.) 
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Figure 1.1.1: 

Let tp : 23# — > R 3 be a sufficiently regular mapping. It is said to be a 
deformation if 

(1.1-2) det(V0)>O 

where is called the deformation gradient and is a matrix given by 



V0 



d<Pi 




dfa 


dX Rl 


dX Rl 


dX R , 


d<p 2 


d(f> 2 


d(f> 2 


dX Ri 


dX R2 


dX Ri 


d(f> 3 


d(f> 3 


dfo 


dX Rl 


dX Rl 


dX Ri 



(pi being the components of (p. 

5 Remark 1.1.1. From d 1.1-21 ) it follows that <p is locally one - one, though 
it may not be globally so. 

The image set 23 = 0(23#) is called the deformed configuration. Note 
that the mapping <p can be written as 

(1.1-3) (p = Id + u 

and the mapping u : —> K 3 is called the displacment. It is also seen 
that 

(1.1-4) V0 = / + Vw 
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where / is the identity matrix and Vw is the displacement gradient. 

The deformation gradient defines the deformation at X = 4>(Xr) up 
to first order. If dt e\ is a line segment parallel to e\ at Xr, it is trans- 
formed into a curve at X whose tangent is dt d\(f>, where d\<$> is the first 
colums vector of V0. The magnitude dt is now 'streched' by dt \d\(p\, 
where |.| stands for the Euclidean norm. The three vectors d\(p, 82$, d^cp 
are independent and, owing to the relation ( I1.1-2L preserve the orienta- 
tion of |<?i,<?2,<?3}- 

It will now be seen how volume, area and line elements are trans- 
formed under the defomation (p. 

(i) Volume elements: The change from a volume element dX^ to dX 
in the deformed configuration comes from the familiar change of 
variable formula in interation theory: 

(1.1-5) dx = det{V<p{X R ))dX R . 

(ii) Surface elements: If dA R is a surface element on &r deformed 
onto a surface elemment dA on 23, then 

(1.1-6) dA = det(V^))|(V<K^r^A R 

where n R is the unit outer normal. (If F is any matrix, F T stands 
for its transpose, F~ l for its inverse and F' T = (F~ l ) T ). 

The formula dl.l-6t will now be proved . This needs some prelim- 
inaries. Let M 3 stand for the set of all 3 x 3 matrices. A tensor will 
be understood simply to be an element of M 3 . 

Let T : 23 — > M 3 be a tensor field. Then its divergence (assuming T 
to be smooth snough) is defined by 

(1.1-7) DIVT = — i-ei. 

dXj 

Thus each component of DIV T is the divergence (in the usual 
sence) of the corresponding row vector of T. By a standard application 
of Green's formula it follows that 

( \ ( \ 

(1.1-8) J DIVTdX= J -Q^ dX e i= J TijnjdA ej = J T n dA 
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where n is the units outer normal to 93. In the same vein DIVr(Tr) on 
tensor fields on 93« can be defined and the analogue of dl.l-Bt can be 
obtained. 

Let T : 23 — > M 3 be a tensor field. ItsPiola Transform is a tensor 
field T R : 93 R -> M 3 given by 

(1.1-9) 7*(X«) - det(V0(X s ))r(X)(V0(X«))- T 

were X = 4>(X R ). 

This is a very useful transformation. The following theorem will 
establish the formula (ll.l-6t . 

7 Theorem 1.1.1. (i) T R {X R )n R dAR = T(X)ndA 

(ii) det(V(/>(XR))(V(l>(XR))- T nRdAR = ndA 

(Hi) det(V0(X R ))|(V0(X R ))- r ««|JA s = dA. 

Proof. It can be shown that (cf. Exercise 1.1-1). □ 



(1.1-10) DIVrT r (Xr) = det(V (f>{X R ))DIVT(X) 

If vr is any arbitrary volume in 23# and # = <p(y R ), then 

J T R {X R )n R dA R = J DIV R T R (X R )dX R 

dv R dv R 

= J det(V ct>(XR))DIVT(X R ))dXR 

Vr 

= J DlVT(X)dX = J T(X)ndA 

Vr dv 

which, as v was arbritrary, proves (i). The assertion (ii) follows by set- 
ting T - I. This is a vector relation and taking the Euclidean norm on 
both sides gives (hi). 
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Remark 1.1.2. The matrix det(V0)(V0) _r = (adj V0)~ r is the matrix of 
cofactors of V0. 

(iii) Line elements: If <f> is smooth enough, cp(X R + SX R ) - (f>(X R ) = 
V<f>(X R )6X R + o(6X R ). 
Thus 

(1.1-11) \<P(X R + SX R )-cj>(X R )\ 2 = 5XlVcf>(X R ) T Vcj>(X R )5X R + o(\5X R \ 2 ) 
which gives the change in length. The matrix 
(1.1-12) C = V(p T V<t> 

is called the (right) Cauchy - Green strain tensor and will play an im- 8 
portant role in the theory. It is used to compute the length of an arc. If 
/(/) is a curve i R in 23#, where / c R is an interval, and I = (p(£ R ) is its 
image in 23, then the length of I is given by 

J Wof)\t)\dt = J ^ Ci j(f(t))f!(t)fj(t)dt 
i i 

where C !; are the components of the matrix C defined above. 

Remark 1.1.3. The matrix 

(1.1-13) B = V<pV<p T 

called the (left) Cauchy-Green strain tensor will be introduced later and 
will play an important role in the constitutive equations. 

Remark 1.1.4. The change in volume depends on a scalar det V0. The 
change is surface elements depends on a matrix, (adj V0) and the change 
in line elements on a matrix, C = V0 r V0. All these will figure in the 
integral representing the energy (cf. Sect. I2.6I ). 

To conclude this section, it will now be examined to what extent the 
strain tensor C is a measure of the deformation. The word 'deformation' 
can be interpreted in two ways - first the formal sense as defined earlier 
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in this section; secondly, in an intuitive way which can be described as 
follows. If <p were merely to consist of a translation and then a rotation 
about a point in space, while it is a deformation in the strict sense, yet 
distances between points are not altered. So intuitively the body has not 
been 'deformed', Such a transformation is called a rigid deformation. 
Thus, (p is said to be a rigid deformation if 

(1.1-14) <f,(X R ) = a + Q(OX R ), 

9 where a e R 3 and Q is an orthogonal matrix whose determinant is +1. 

The vector a above represents a translation and the matrix Q a rota- 
tion. The following notation will used for various classes of matrices: 

M 3 = j/c idel(/- ) > 0) 
Q 3 = {F e M 3 \F T F = FF T = 1} 
O 3 = {FeG 3 |det(F) = +l} 

•:' - (/- c !•: i/' y = f} 

S 3 = {F e S 3 |Fis positive definite}. 

Thus Q e O 3 . Observe that if cp is rigid then C = Q T Q = L In fact, 
under suitable hypotheses, the converse is also true. 

Theorem 1.1.2. Let Q. be an open connected subset of R 3 . Let <p e 
C l (£l\ W?) such that for all x e Q, 

(1.1-15) V<p(x) T V(p{x) - / 

Then, there exists a vector a e R 3 and a matrix Q € O 3 such that, 
for all x € O 

(1.1-16) <f>(x) - a + Q(0x). 

Proof. Cf. Exercise 1.1-2 □ 

Theorem 1.1.3. Let O be an open connected subset ofR 3 and let <p,\p € 
C 1 ^; R 3 ) such that for all x e Q. 



(1.1-17) 



V(p{x) T Vil/{x) - V<Kx) r V0(x). 
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Assume further that ip is one - one and that det(Vt/r(x)) + Ofor all 
jefl. Then there exists a e K 3 and Q e O 3 such that for all x e Q. 



Proof. Consider the mapping 9 - <po\p on ^(Q). Clearly is con- 
nected. Also, under the given conditions, it is open by the theorem of 
invariance of domain. Further, 6 e C l (ifr(Q.); R 3 ). Now from (11.1-171 it 
follows that 6 satisfies dl.l-15t and so the previous theorem applies to 8 



Thus if two deformations have the same strain tensor then, upto a 
rigid deformation, they are the same. Thus C 'measures' the 'deforma- 
tion' upto a rigid transformation. Naturally, a measure of the deviation 
from a rigid deformation is obtained from C — I. The Green-St Venant 
strain tensor, E, is defined by the relation 



(1.1-18) 



(p{x) — a + Q\jj{x). 



and the result follows. 



□ 



(1.1-19) 



C-I = 2E 



In terms of the displacement gradient, 



/ + 2E = C = V(p T V(f) = I + Vu T + Vu + Vu T Vu 



or, componentwise, 



(1.1-20) 



Ey = -(djUj + djUi + diU m djU m ) 



where 5; stands for 



d 



dX Ri 



Exercises 



1.1-1 Prove the Piola identity 



DIV R (det{Vcf>{X R ))(V4>({X R )y T ) = 0. 



Deduce the relation dl.l-lOt from this. 
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1.1-2. Prove Theorem 1 1 . 1 .21 (Hint: First show that at least locally, (p is 
an isometry; then show V0 is locally contant and use the con- 
nectendness of CI.) 

1.1-3. Let <p : W — > W,n > 2, be continuous. Assume that there exists 
€ > such that for all x,y e R" with \x-y\= {, \<p(x) - <p(y)\ = t. 
Show that (p is an isometry, i.e. there exists a € R" and Q e O" 
such that for all x € R.' 1 



1.1-4. Given a tensor field r : Q. — > S 3 , find necessary and sufficient 
conditions such that there exists a mapping cf> : Q — > R" with 



1.2 Euilibrium Equations 

The equilibrium equations give the relationship between the given forces 
acting on a body and the state of "stress" (to be defined below) which 
results as a consequence of these forces. 

Let the mass density as X e 23 be given by p(X) while that at Xr e 
23# is given by pr{Xr). The applied forces in 23 are of two kinds. 



0(x) = a + Qx. 



r = V(f> T V<f> 




h{X)dA 



n 



Figure 1.2.1: 
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(i) Body for volumic) forces: b : 93 — > K 3 . The elementary force on 
a volume element dX will thus be p(X)b(X)dX. An example of a 
body force is gravity and in this case b = (o, o, -g). 

(ii) Applied surface forces: t\ : 323i — > K 3 , where 523i is a portion of 
the boundary <993. If ciA is a surface element, the applied force on 
it will be t\dA. An example of a surface force is a pressure load 
where t\ = -pn, p € R, n the normal to dA. 

A system of forces in 93 consists of body forces (identical to (i) 
above) and surface forces t : 23 x Si — > K 3 where Si is the unit sphere 
in K 3 , i.e., 

Si = {xeR 3 ||x| = 1}. 

If ■& is any sub volume of 93, dA a surface element of d& and n the 
normal to it, the surface force t(X, n)dA acts in it. Note that this is inde- 
pendent of &, i.e. if #i were another subvolume and dA lay on dv\ with 
the same n as normal, the force acting on it will remain as t(X, n)dA. 13 
Further if dA c <993 and n were also normal to 523, it is required that 

(1.2-1) t(X,n) = h(X). 

The vector t(X, n) is called the Cauchy stress vector. 
The following axiom is the basis of Continuum Mechanics in gen- 
eral, and consequently of the theory of elasticity in particular. 

AXIOM OF STATIC EQUILIBRIUM. Let 13 be a deformed configu- 
ration in static equilibrium. There exists a system of forces such that for 
any subdomain & c 23,f/ie corresponding system of forces is equivalent 
to zero (in the sense of torsors). Thus 

(1.2-2) J p(X)b(X)dX + J t(X,n)dA = o. 

(1.2-3) J OXAp(X)b(X)dX + J OXKt(x,n)dA = o. 



The wedge A stands for the usual cross product of vectors in K 3 . 
The following notation will be useful in manipulating cross products. 
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For indices i, j, k taking values 1, 2, 3 the tensor of rank 3, e^, is 
denned by 



(1.2-4) e ijk 



+ 1 if(i, j, k)is an even permutation of(l, 2, 3), 
-1 if it is an odd permutation of(l, 2, 3), 
otherwise. 



Then for vector a, beM? , 

(1.2-5) ahb = eijkdjbkei. 

14 The following consequence of the axiom of static equilibrium is of 

paramount importance. 

Theorem 1.2.1 (Cauchy's Theorem). Let pfC°(23 ; R), &eC°(23;R 3 ), 
?(., ri) eC x (^&\ R 3 ) and t(X, .)eC°Q] i ; R 3 )- Then there exists a tensor field 
TeC 1 (23; M 3 ) such that 

(1.2-6) t(X, n) = T(X)n, for all Xe23, neSi , 

(1.2-7) DIVT{X) + p(X) b(X) - 0, for all Xe93, 
(1.2-8) T(X) - T T (X), for all Xe93. 

Proof. Let Xq be any point in 23. Consider a tetrahedron ■& with vertices 
Xa. Vi . Vi. Vi. as shown in Fig. 11.2.21 □ 




Let n^eYji be the normal to the plane V\ V2V3 and keep n fixed to 
begin with. Let the distance of X to the plane be 8, let Sj{i - 1, 2, 3) be 
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the surface opposite the vertex Vi(i - 1, 2, 3) and S the surface opposite 
Xq. Since p, b are continuous on 58, they are bounded. Thus by dl.2.21) . 

| J t(x,n)dA\ < KVol{&) 
ao 

K being a constant independent of 8. Since Vol(#) = K\6 3 , A(5) = Area 15 
of S = Kid , Ki,K2 being independent of 8, it follows that 

(1.2-9) lim — \ t{X,n)dA = 

9& 



Now, 



— — r ?(x, f 



lim — — n) dA = t(X , n ) 

s->o A{5) 

s 



and 



hmS -> 0— - f *(X,n)dA - (n .ei)t(X , -eft 
A(o) J 

s, 

using the continuity of the given functions. Hence by ( 11.2-91 . 
(1.2-10) r(X ,n ) - -(n .ei)t(X , -e ). 

If ti — > e/, again by continuity of t it follows that 
t{X ,ej) = -t(X ,-ej). 

Thus, on substituting this in dl.2-10l >. 
(1.2-11) t(X ,no) = t(X ,ej)nj. 

Setting 

f(X ,e y ) = Tij(X )ei 

the equation (ll.2-6t follows. The smoothness of T results from that of t 
w.r.t. X. 
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Using dl.2-61 ) in ( 11.2-21 ). for any volume &, 

= J p(X)b(X)dX + J T(X)ndA 

= J p(X)b(X) + DIV(T)dX, 



16 from which (11.2-71) follows as # was aribitrary. 
Finally by and (fl~2~3t 

= J ey jfc X; P (X)fc jfc (X)^ + J 



T k( n e dA 



f 5X / r 



dX+ | eij k XjT kf ncdA 



iCkTke, 



using dl.2-7t . Since # was arbitrary, 
which is just a restatement of d 1 .2-8b - 

Remark 1.2.1. Given a tensor field T : 93 — > M 3 satisfying (11.2-71) and 
the vector field t(X, n) = T(X)n satisfies SFT3k and (IT331) . 

The tensor T(X) obtained in the above theorem is called the Cauchy 
stress tensor at the point Xe23. 



e 3 

4 



Tne 



-C2 



ei 




e 1 )=T il e i 



Figure 1.2.3: 
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Remark 1.2.2. The components of T can be interpreted as follows. If 17 
an element dA has normal e\ then the Cauchy stress vector acting on it, 
t(X, e\) has components T\\, Ti\ and T$i and so on. 

The Cauchy stress tensor thus satisfies a boundary value problem: 

DlVT+pb = Q 



T = T T \ 



in SB 
T n - t\ on 593i 



Let u.v stand for the usual scalar product in R 3 , i.e. v.v = m,-v,-. If 
A, BeM 3 , denote 

(1.2-12) A : B - AijBij = tr(AB r ). 

This is an inner product in M 3 with the associated norm 

(1.2-13) ||A|| = Va~A~ 

Using Green's formula, a variational form of the boundary value 
problem can be obtained. 

If T is a tensor field and <Qj is a vector field on 25 then 

J DIV T.@dX = J ^®idX 

SB 33 

= ~I Tij 7m dx + J Ti J® in J dA 

SB dSB 

= - Jt : GRAD@dX + J" 



G/?AD@dX + J r«.0dA. 

SB 3 SB 

In particular, if T is a solution of the above boundary value problem 18 
and if (Q) vanishes on <923o = <923\<993i, then Green's formula above gives 



f 



- (DIVT + pb).®dX 
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= J(-T : GRAD® + pb.®)dX + J h.®dA. 

SB <9SB] 

Conversely if the above relation is satisfied for all H) vanishing on 
(923o then T is a solution of the boundary value problem. Thus 

Theorem 1.2.2. The following are equivalent: 



(i) 



(DIVT +pb = 0, inSB 
I Tn = t\ in 323i 



(ii) For all <H> : 93 -> R 3 , © vanishing on <923 , 

(1.2-14) J~ r:G/?AD@t/X- J" pb.®dX + J h.®dA. 



The equations ( 11.2-141 form the so-called variantional formulation 
of the boundary value problelm (i). In Mechanics, it is also known as 
the Principle of Virtual Work in the deformed configuration. 

The equations of equilibrium were established in the Eulerian vari- 
able, X, in the deformed configuration. However, this is of no use for 
computation as the deformation <p is unknown. So, the equations must 
be written in the reference configuration, which is a fixed domain given 
a priori, in terms of the Lagrangian variable, Xr. In doing this, it is 
desirable to retain as much of the divergence form of the equations as 
possible so that a similar variational formulation can be obtained in the 
reference configuration. It is here that the merit of the Piola transform 
is seem. 

19 The Piola transform of the Cauchy stress tensor T, called the. first 

Piola-Kirchhoff stress tensor, is denoted by Tr. Thus surya 



T R = det(Vcf>)T{VcPT T ■ 
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Figure 1.2.4: 

By the principle of conservation of mass, it is known that 
PR (X R )dX R = p(X)dX. 

By defining 

(1.2-15) b R (X R ) = b{X),orb R = bo<p 

it follows that 

p R b R dX R = pbdX. 
Note that, a priori, b R depends on (p. 

Remark 1.2.3. Since it is known that dX - det(V cf>)dX R , it follows that 

(1.2-16) p(X) = -^g-T- 

det \(p(X R ) 

Since the density at any point (in either configuration) has to be 20 
finite and positive, this, if not any other, is a necessary reason for a 
deformation to satisfy. det(V^) + 0. 

Multiplying equation dl.2-7t by det (V0) on both sides, it follows 

that 

(1.2-17) DIV R T R + p R b R = in 23# 

Thus the divergence form is preserved. Note however that T R is not 
symemetric. A symmetric tensor to T R can be defined. It is the second 
Piola-Kirchhoff stress tensor, Yj R , given by 

(1.2-18) Yj = det(V0)(V0)- 1 r(V0)- r 
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It is related to T R by 
(1.2-19) J] = W)" 1 ^. 

R 

Remark 1.2.4. It is understandable that T R is not symmetric as it be- 
longs partly to the reference configuration and partly to the deformed 
configuration and symmetry does not make much sense in such a situa- 
tion. 

Now we turn to the transformation of the surface forces. Theirs? 
Piola-Kirchhoff stress vector is defined so that 

(1.2-20) t R (X R ,n R ) = T R (X R )n R . 




n 



Figure 1.2.5: 

21 Recall that T R (X R )n R dA R = T(X)ndA and so 

t R (X R ,n R )dA R = t(X,n)dA. 

If d^B\ R is the portion of <923# mapped by <p onto <923i, define t\ R : 
(993 \ R — > K 3 by t\ R dA R = t\dA. Again, a priori, fm depends on <p. 
Explicitly, by Theorem ll.l.ll 

(1.2-21) t mXR) = det(V0(X R ))|(V0(X s ))- r « R k 1 (0(X R )). 

The following result is easy to establish. 

Theorem 1.2.3. The equilibrium equations in the reference configura- 
tion are given by 

( 1 .2-22) DIV R T R + p R b R = oin% R 



1 .2. Euilibrium Equations 
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(1.2-23) (Vcf>)Tl = T R (V<f>f in^ R 

(1.2-24) T R n R = t R ond% R . 

Equivcilently, in terms ofYjR 

(1.2-25) DIV R (V<p J) + PRb R = in » R 

R 

T 

(1.2-26) Yj = Yj in * R 

R R 

(1.2-27) V<pY j n R = t m ond% m . 

R 

Again, this is equivalent to the variational equations 
(1.2-28) jT R :V6dX R = j p R b R .6dX R + J t lR .GdA R 

/or a// : S8« — > K 3 vanishing on d%i oR = <323r\<923ir. 

Remark 1.2.5. Equations (ll.2-28t go under the name of the principle of 22 
virtual work in the reference configuration. 

To conclude this section, some classes of applied forces are consid- 
ered. Recall that while p R is completely known, b R and ?ir depend in 
general on <p which is unknown. 

A body force (resp. applied surfaces force) is a dead load if b R (resp. 
t\ R ) is a function of X R only, independent of (p. 

An example of a body force which is a dead load is gravity which 
is constant; b = (o, o, -g). A trivial example of an applied surface force 
which is a dead load is t\ =0! The pressure is an example of an applied 
surface force which is not a dead load: 

(1.2-29) t x = ~pn 

where p > o indicates an inward directed force (pressure) and p < 
indicates one which is directed outward (traction). Now 

t lR = -pdet(V(p)(V(l>y T n R on <923 R 
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which clearly depends on </>\ 

A body force is said to be conservative if there exists a function 
{3 : R 3 x 23 R -> {3(<p,X R )eR, such that 

(1.2-30) b R (X R ) = VfPMXn), X R ), 

for all X#e23« and all deformations (p. If which is the case then 

(1.2-31) J" p R b R .6dX R = B{<p){6) 

where 

(1.2-32) fi(<A) = J PR (X R )j3ty(X R ),X R )dX R . 

SB* 

23 A body force which is a deal load is conservative, fi(Q,X R ) = b R (X R ). 
O. 

An applied surface force is conservative if there exists a function 
Ti : R 3 x <923i R -> Ti(<f>,X R )eR such that 

(1-2-33) ?i R (X«) = V$Ti(0(X R ),X R ). 

Then again 

(1.2-34) J t 1R .6dA R = T[(<f>)(6) 



where 



(1.2-35) W)= J Ti(^ R ),X R )JA R . 

Ara applied surface force which is a deal load is conservative; 
X\(<p,X R ) = tm(X R ).0. A pressure load is conservative (Exercise 1.2-3). 



Exercises 



1.3. Constitutive Equations 
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1.2-1 . (Da Silva's Theorem). Given any system of applied forces (with 
323 1 = 323) show that there exists QeO 3 such that 

SB as 

J p(X)Q T (OX)Ab(X)dX + J Q T (OX)At(X)dA = o. 



as 



How many solutions exist? 

1.2-2. Show that the fundamental axiom of static equilibrium is equiv- 
alent to 

J p(X)b(X).v(X)dX + J t(X,n).v{X)dX = 

for every volume ft c 23 and for every infinitesimal rigid dis- 24 
placement v, i.e., 

v{x) = a + bAOx, a, beR 3 . 

This is sometimes also called the principle of virtual work. 1.2-3. 
1.2-3 Show that a pressure load is conservative. 

1.3 Constitutive Equations 

Given a body acted on by a system of forces, one's main objective is 
to compute the deformation (f> which has 3 component functions. As 
a natural intermediary, the stress tensor T has come in which has 6 
components (taking into account its symmetry). But so far, the boundary 
value problem obtained via the equilibrium equations has yielded only 
3 equations (cf. (I1.2-7I) ). Thus 6 more equations must be found. 

From the physical point of view, observe that in obtaining the equi- 
librium equations, no property of the material under consideration has 
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been used. Since different materials react differently to the same forces, 
obviously these equations alone cannot describe the response of the ma- 
terial. 

Thus one is led to finding more equations to complete the system. A 
material is said to be elastic if there exists a mapping 

f : FeM 3 + -> f (F)eS 3 

25 such that for any deformed configuration and any point X = <p(X R ), 

(1.3-1) T(X) = f(VO(X«)). 

The map f is called the response function and (11.3-11) is called a 

constitutive equation. 

Remark 1.3.1. The map f above does not depends explicitly on X R . 
Such that a material is called homogeneous. If it were that 

T(X) = f(X R ,V®(X R )) 

the material would be called a non-homogeneous elastic material. 

If T R is the Piola transform of T then it follows that 

(1.3-2) T R = det(V0)f (V0XV0)- 1 d = f R (Vc/>) 

which gives a reponse fucntion f R : — > M 3 for T R . Similarly it is 
possible to write one for ~L R in terms of a response function E# : M 3 + — > 

S 3 . 

Theorem 1.3.1 (Polar Factorisation). Let F be an invertible matrix. 
Then there exist an orthogonal matrix R and symmetric, positive defi- 
nite matrices U and V such that 



(1.3-3) F = RU = 

Such a factorization is unique. 
Proof. Cf. Exircise 1.3-1. 



VR. 



□ 



1.3. Constitutive Equations 
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Remark 1.3.1'. If F € M 3 then R 6, Q 3 . If G € § 3 there exists a unique 
matrix H e S> such that // 2 = G. It is usual to write H = G 1/2 . It can be 
seen that U = (F T F) 1 ^ 2 and V = (FF T ) 1 ^ 2 , in the above theorem. Since 
V = RUR T , U and V are similar. Then so are B - FF T and C = F T F. 

The constitutive equation ( 11.3-U can be written componentwise as 



ever, the functions f, y - cannot be chosen arbitrarily. They must somehow 
reflect an intrinsic property of the material in equation, irrespective of 
the coordinate system chosen. This is the idea embodying the 

AXIOM OF MATERIAL FRAME INDIFFERENCE. The Cauchy stress 
vector t(X, n) = T(X)n should be independent of the particular basis in 
which the constitutive equation is expressed. 

Theorem 1.3.2. The following are equivalent. 

(i) A response function t : M+ — » S 3 satisfies the axiom of material 
frame indifference. 

(ii) For every Q € O 3 and for every F € M 3 , 




(1.3-4) 



t{QF) = Qt{F)Q T . 



(iii) For every F e M 3 



if F — RU is its polar factorisation then 



(1.3-5) 



f(F) = RT(U)R 



,T 



(iv) There exists a map 2,r : — > S 3 such that 



(1.3-6) 



1 R (F) = 1 R {F T F) 



for every F € M 3 



Proof, (i) o (ii) Instead of rotating the coordinate axes the same effect 
can be achived by rotating the deformed configuration. □ 
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Figure 1.3.1: 

Rotating S by a map Q e Q\, let X map into X' . The normal n at 
any goes to Qn and t{X, n) goes to Qt(X, n). Thus 

t(X', Qn) = T'{X')Qn 
t{X',Qn) = Qt{X,n) = QT(X)n. 



28 Since n is arbitrary, it follows that 

(1.3-7) T'{X') = QT{x)Q T 

Thus 

f (QVcf>(X R )) = Qf(V<f>(X R ))Q T 

for any Q e 3 + and any F = V(f> e M 3 + . This shown that (i) => (ii) 
Simply retracting the argument proves the converse, 
(ii) o (iii). If F = RU, then by (ii), since R e 3 + 

f{RU) = RT(U)R T 

which is (iii). Conversely, assuming (iii), if F = RU then the polar 
factorization QF is (QR)U for Q e O^, as the factorization is unique. 
Thus 

f(QF) = QRf(U)R T Q T = Qf(F)Q T . 



1.3. Constitutive Equations 
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(ii) o (vi). Since F - RU implies U = (F T F) l/2 , 

f(F) = Rf(U)R T 

= FU^t{U)U' l F T 
= FS (F T F)F T 



where S : S> — > S 3 . Conversely, if f(F) is of the above form, then if 
F - RU, 

f(U) = US(U 2 )U 

and 29 

f(F) - FS (F T F)F T 
= FS(U 2 )F T 
= FU~ l f{U)U' l F T 
- RT(U)R T . 

Now 

S R (F) = det(F)F" 1 f(F)F" :r 

= (det(F T F)) 1/2 5(F T F) - 1 R (F T F). 

Remark 1.3.2. If one of the response functions, say T, can be written of 
either variables F, F T F = C, FF T = B or E ( where C - I + 2E), the 
following notation will be employed when the different dependences are 
expressed: 

T - f(F) = T(F T F) - T(FF T ) - P(£) 

In the above theorem it has been proved that it is enough to know 
the action of T on a relatively small class of matrices like §>. 

A material or response function is said to be isotropic if the Cauchy 
strees tensor (or vector) computed at a given point in the deformed con- 
figuration is the same if the same if the reference configuration is rotated 
by any rigid defomation. 
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While the axiom of material frame indifference is an axiom to ver- 
ified by any response fucntion, isotropy is a property of a particular 
material. There can be materials which are non-isotropic; for instance, 
a body up of layers of different materials. 

30 Theorem 1.3.3. The following are equivalent. 

(i) A response function f : M 3 — > § 3 is isotropic. 

(ii) For every F € M 3 and for every Q € G 3 , 
(1.3-8) f(F) = f(FQ). 

(iii) There exists a map f : S 3 — > S 3 such that for every F € M 3 , 
(1.3-9) f(F) = T(FF T ). 

Proof, (i) o (ii). Let Q € O 3 . Rotate the reference configuration about 
a point Xr so that if Xr € then 

(□) 0(X R ) = X r + Q t (X r X r ). 

Then 




Figure 1.3.2: 

The response function is isotropic if and only if 

f (X) = f(V<f>(X R )) = f(Vf(X R )) 



1.3. Constitutive Equations 
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i.e., 



f (V<p(X R )) = f (V<p(X R )Q). 



31 (ii) o (iii). Let FF T = GG T ,F,G e M 3 . Then G l F e 3 . Hence by 



So it is clear that f(F) depends only on FF T . Conversely if f(F) - 
T(FF T ) then for Q e 3 + , 



Remark 1.3.3. By the axiom of material frame indifference, the consti- 
tutive equation could be expreesed in terms of a function of C = F T F 
and this involved rotating the <fe/bn«e<iconfiguration 23. By isotropy, 
the same could be expressed in terms of a function of B = FF T and 
this involved rotating the reference configuration 23#. Thus these two 
nations seem to be 'dual' ot each other. 

Remark 1.3.4. For non-isotropic materials it can be shown that 

f{F) = f{FQ) 
for all F e but Q variying over a subgroup of 0+. 

In what follows, the material will allways be assumed to be isotropic. 

Before proving a very powerful and elgent result on the structure of 
a reponse function which is isotrophic and material frame-indifferent, 
the following definition is needed. 

Let A e M 3 . Define i A to be the triple (n(A), 12(A), 13(A)) where 
11(A), 12(A) and 13(A) are the principal invariants of A, and 



(ii) 



f(G) - f(G(G~ l F)) = f(F). 



f(FQ) = T(FQQ T F T ) - T(FF T ) - f(F). 



□ 



(1.3-10) 



det(A - AI) - -A 3 + n(A)A 2 - i 2 (A)A + 13(A). 



If A = (aij) and Ai,A 2 , A3 are its eigenvalues, then 



(1.3-11) 



i\(A) - an - tr(A) - A\ + A 2 + A3 . 
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12(A) = -(a iiaj j - ajjOij) = ~((tr(A)) 2 - tr(A 2 )) 

(1.3-12) = tr(adj A) = A X A 2 + A 2 A 3 + A 3 A\. 

(1.3-13) 

13(A) = det(A) = \((tr(A)f - 3tr(A)tr(A 2 ) + 2tr(A 3 )) = AiA 2 A 3 . 
6 

Further, if A is invertible, 

(1.3-14) 12(A) - (det A)tr(A~ l ). 

The following theorem is one of the most important results in the 
theorey of elasticity. 

Theorem 1.3.4 (Rivlin-Ericksen Theorem). A response function t : 
M? + — » S is isotropic and material frame indifferent if and only if 
it is of the f(F) - f(FF T ) where the mapping f : S> — > S 3 is of the 
form 

(1.3-15) f(B) = p (i B )I + ^(i B )B +/3 2 (ib)B 2 

for all B € §>. where B ,B\,B 2 are real valued functions. 

Proof, (i) Let f : M 3 + — > S 3 be material frame indifferent and isotropic. 
Then by isotropy f(F) - f(FF T ) for some mapping T : S 3 — > S 3 . Let 
Q e O 3 and B e S 3 . On one hand , by isotrophy 

f(QB l/2 ) = T(QB l/2 B 1/2 Q T ) = T(QBQ T ). 

On the other hand, by the material frame indifference, 

f(QB l/2 ) = Qf(B l/2 )Q T 

= QT(B l/2 B l/2 )Q T = QT(B)Q. 



Thus T satisfies , for all Q e O 3 , and B e S 3 
(1.3-16) T(QBQ T ) - QT(B)Q T . 



1.3. Constitutive Equations 
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Conversly, let r : § 3 — > § 3 satisfy 11.3-161 and let f(F) = T(F T F). 
Then clearly, T is isotrophic. If Q e O 3 , then 

f(QF) = T(QFF T Q T ) = QT(FF T )QT. 
= Qf{F)Q T 

and so 7* is material frame indefferent. 

Thus it is now enough to check that a mapping t : § 3 — > § 3 sat- 
isfying (ll.3-16t is of the form d 1 .3- 151) - (The converse is immediate to 
varify). 

(ii) Let T : § 3 -> § 3 varify 11.3-161 It will now be shown that 
any matrix which diagonalizes B € S 3 also diagnalizes f (B), i.e., any 
eigenvector of B is an eigenvector of T (B). 

Let B € S 3 and Q € O 3 (we can always assume that) such that 

Q T BQ = diag (^) 

where X\,X2, ^3 are the eigenvalue of B. Define 

£1 = diag (1, -1,-1), & = diag (-1, 1, -1), Q 3 = diag (-1, -1, 1). 

Then Q k e0 3 + ,k= 1,2,3. 
Also, 



So, 



Q T k Q T BQQ k = diag Xi = Q T BQ. 



QlQ T f(B)QQ k = f(Q T k Q T BQQ k ) 



= f(Q T BQ) 
= Q T T(B)Q. 



If D - Q T T(B)Q, then 



Q T k DQ k = D,k = 1,2,3- 

If the diagonal entries of Q k are q k (= 1 if i = k, -1 if i + k), then it 
follows that 

D u = q^Dijq k j for all 1 < i, j, k < 3. 
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Thus if i - k + j, then 

D kj = -D kj or D kj = o. 

Hence D is diagonal and this proves the claim. 

(iii) It will now be a shown that if f satisfies dl.3-16t then, for all 



Be 



:3 



(1.3-17) T(B) = b Q {B)I + b { {B)B + b 2 {B)B 2 , 



b a , a = 0, 1, 2 being real valued functions on 



;3 



35 Case 1. B has 3 distanct eigenvalues A\, ^2,^3 with corresponding or- 
thonormal eigenvectors pi, P2, P3- Then 

(1.3-18) I = P\p\ + PipI +P3PI 

(1.3-19) B = A yPl p[ + A 2 P2P T 2 +A 3 p 3 pl 

(1.3-20) B 2 - A\ PlP \ + A&zpl + Ajpipl 

Since the A, are distinct, the Vandermonde determinant 



det 



1 1 1 

A\ A 2 A?, 
A\ A\ A\ 



is non-zero and so in S 3 , the span of Ptpf , i — 1, 2, 3 is equal to that of 
I,B,B 2 . But f(B), by step (ii) above, has the same eigenvectors as B. 
So 

(1.3-21) f(B) = p.\p\p[ + p 2 p 2 p T 2 + P3P3PI 

which implies that T(B) € span {/, B, B 2 }. 

Case 2. A\ + A 2 - A 3 . Again one can write dl.3-18t and dl.3-191 ). Then 
the span of p\p\ and p 2 p\ + P3p\ is that of / and B. By step (ii), it can 
be seen that \l% — P3, since any non-zero vector spanned by p 2 and p 3 is 
also an eigenvector for T (B). Thus in A1.3-2H 

f(B) = p.\p\p[ + p. 2 {p 2 p\ + P3p\) 
which shows T (B) € span (/, B). 



1.3. Constitutive Equations 
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Case 3. A\ — A% — A3. In this case, one can similarly see that B, T(B) 
are both scalar multiples of /. 

(iv) Case. [TJ A\,A 2 ,A^ are distinct eigenvalues of B € S>, . 
Let Q e Q 3 . 

f(QBQ T ) = b (QBQ T )I + bi(QBQ T )QBQ T + b 2 (QBQ T )QB 2 Q T 
= Q(b (QBQ T )I + b\{QBQ T )B + b 2 (QBQ T )B 2 )Q T 

But 

f(QBQ T ) = QT(B)Q T = Q(b () (B)I + b,{B)B + b 2 (B)B 2 )Q T . 
Thus b a : §> —* R, a = 0, 1 , 2 satisfy the functional identity 
(1-3-22) b a (QBQ T ) = b a (B) 

for all B € § 3 and for all Q € O+. Thus if Q diagnalizes B, it is seen 
that such a function b a must be a function of the eigenvalues of B only. 
Now choosing Qi e O+, z = 1, 2, 3 as 








1 







-1 
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Gi = 


1 








,Qi = 








1 


,Q3 = 





-1 













-1 







1 







1 









it is seen from H3-221 that b a is a symmetric function of Ai,A 2 , A3, i.e,. 
b a (B) = B a {t B ). 

This proves the theorem completely. 

Theorem 1.3.5. (a) Given 23# and an isotropic material frame indif- 
ferent material, then in any deformed configuration 23 - 0(23«), the 
Cauchy stress tensor is given by 

(1.3-23) T(X) = f (V0(X*)) = f (V0(X R )V<£(X«) r ), 

f : Si-* S 3 satisfying (I1.3-15H . 

(fr) 77ze second Piola-Kirchhoff stress tensor is given by 
(1.3-24) S R (X S ) - £ R (y<KX R )) = I«(V0(X R ) r V0(X R )) 
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where Er : Sf. — » S 3 satisfies, for all C € § 3 , 

(1.3-25) f R (C) = y {i c )I + 7l (i c )C + y 2 (i c )C 2 . 

Proof. Observe that 

£ R (F) = det(F)F~ l f(F)F- T 

= {det(F T F)) l/2 F~ l T{FF T )F~ T 

= (detQ^F-^oOg)/ + J 8i(wr B )5 + B 2 {i B )B 2 ]F- T 

where C = F T F, B = FF T . But these are similar. So ig = i c . 
Further 

F- l F~ T = C _1 

F l BF T = / 
F _i g 2 F -T _ c 

By the Cayley-Hamilton theorem, 

-C 3 + ri(C)C 2 - i 2 (C)C + j 3 (C)7 - 
or C" 1 = -^-(C 2 - n(C)C + 12(C)/) 

38 where 23(C) = det C ^ 0. Thus it is clear from these considerations that 
Sj? can be expressed in terms of C as in ( 1 1.3-241 ) - dl.3-25t . □ 

It was seen in Section 11.11 that the Green-St Venant strain tensor 
E, given by C = / + 2E, 'measures' the actual deformation. If E# 
is sufficiently smooth it is possible to express it in terms of E. More 
precisely, the following result is true. 

Theorem 1.3.6. Let 23# be the reference configuration of an isotropic, 
material frame indifferent elastic material. Assume that the functions 
y a , a — 0,1,2 of ( 11.3-251 ) are differentiable at ij - (3, 3, 1). Then 

(1.3-26) T. R = £ R (I + IE) = -pi + (A{trE)I + 2pE) + 0(E), 

where p, A and p are constants. 



1.3. Constitutive Equations 
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Proof. Using the relations 

tr(C) - 3 + 2tr(£) 
tr(C 2 ) - 3 + 4tr(£) + o(E) 
tt(C 3 ) = 3 + 6tr(E) + o(E) 

and the relations (11.3-llb - dl.3-13b . it follows that 

H(Q = 3 + 2tr(£) 

i 2 (Q = 3 + 4tr(£) + o(E) 

i 3 (C) = l + 2tr(E) + o(E) 

so that 

J{ic) = Tin) + (2^0/) + 4^(z/) + 2^0/)) tr(E) + 0(E) 

Oil Oil o«3 

where y = (70,71,72)- This yields dl.3-26t . In particular 39 
(1.3-27) p = -( 7l (u) + jy(ly) + 7200). 

(1.3-28) ^ = 2(2^(1!) +4^ + 2^)). 

(1.3-29) /i = 71O1) + 272O1). 

□ 

A reference configuration is a natural state if 'there is no stress in 
it', i.e., p = 0. In this case 

(1.3-30) Lr = £j(£) = itr(£)/ + 2pE + o(E) 

and /I and p are called Lame's constants. It is possible to obtain a priori 
some information on the nature of the Lame's constants. 

Let 23# be a natural state and have a 'simple form'. Let cf) e : 23# — > 
R 3 be of the form 

(1.3-31) <f (X R ) = X R + € h(Xk) + o(e;X R ), 
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where e> o is a small parameter, and where Vu{Xr) = G, a constant ma- 
trix. Such a deformation is a special case of the so-called homogeneous 
deformations where V0 is a constant vector. Then 

T e (X) = f(I+ € G + o(e;X)) 

= f(I+ € G) + 0(e;X),X = 0(X R ). 

Thus DIVT e (X) = o(e;X) i. e. to within first order in e, there can 
be no body force: such deformations can only be produced by applied 
surface forces. Now it can be seen that 
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(1.3-32) 



T e {X) =g {A(tr G)I + viG 1 + G)) + o(e; X) 



using the fact that 23^ is a natural state. For particular <p e considered 
the corresponding T e is of some simple form. This can be substituted in 
( 11.3-321 ) and it is thus possible to obtain inequalities for A and fi. 

Experiment 1. Let Br be a rectangular block. Choose 



(1.3-33) 



clef 



u\X R ) =eu(X R ) =g 



Xr 2 





Thus the body deforms as shown in figure [1331 

T 12 e l 




Figure 1.3.3: 
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Then it is logical to assume T% =€ Tn + o(€) where Tu > o. This 
follows from the interpretation of the components of the stress tensor 
(cf. Remark fl.2.2t . Comparing this form with (ll.3-32t if follows that 



(1.3-34) 



li > 0. 
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Experiment 2. Let 23# be a sphere which is contracted by means of a 
normal pressure. Thus 



(1.3-35) 



u%X R ) =6 



+ o(e;X R ). 



Thus 
(1.3-36) 



T e (X) = -peI + o(e; X), p > 0. 




Figure 1.3.4: 



It can then be shown that 
(1.3-37) -/? e / = - e (3A + 2fx)I + 0(e) 
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from which it follows that 

(1.3-38) 3A + 2p>0. 

Remark 1.3.5. This precludes incompressible materials! An example 
of an incompressible material is rubber. 



Experiment 3. Let 23# be a cylinder which is stretched as in figure [l~3~5l 



d = (1 — ve)d,R 



't E (X,d 3 ) = 







<X~^) 
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XR "el 











• s 










liR 


Pr~ 


P - 














































Fi 


sure 


1.3.5: 



h = (l + e)h R 
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Now 



(1.3-39) 



and 



(1.3-40) 



u%X R ) =e 



-vXri 

-VX R1 
Xr3 



+ 0(e;X R ),v> 



T e (X) =e 





£ 



+ 0(e;X). 
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It can now be shown that 

A u(3A + 2u) 

(1.3-41) v= ,E=— — 

2(A + n) A+n 

Since fj. > and 3A + 2/i> 0, it follows that A + fi > 0. Since v > 
it follows that 

(1.3-42) A > 0. 

Thus /t > o and /u > 0. (This does not make sense for incompress- 
ible materials). The number v is known as Poisson's ratio and E as 
Young's modulus. The Lame's constants can be expressed in terms of 
these quantities: 

Ev E 

(1.3-43) A = 



(l + v)(l-2vr 2(1 +v) 
Thus A > O and fi > O is equivalent to 

(1.3-44) 0<v<l/2,£>0. 

(For an incompressible material, v = 1/2). 

An elastic material is said to be a St Venant-Kirchhoff material if 44 

(1.3-45) Z* R (E) = A{trE)I + 2fiE. 

It is also expressible in terms of C: 

(1.3-46) f R (C) = &n(C) - 3)-fAl + nC 



.2 

Then the Cauchy stress tensor can be written as 

A 



(1.3-47) T = T(B) = fe(fi)) 1/2 {-(*i(fi) - 3)-^B + ti(n(B)) i/2 B 2 . 

Thus such a material is isotropic and material frame indifferent (cf. 
Theorem ll.3.4t . 
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1. Description of Three - Dimensional Elasticity 



Remark 1.3.6. While the relation dl.3-45t between Er and E is linear, 
as a function of u, Er is non-linear since the dependence of E on u non- 
linear (cf. (I1.1-20H ). 

The relation dl.3-45t can be written componentwise as follows: 

Ei? y = AE kk 6ij + 2/uEij 
(1-3-48) d = a ijk[ E ke 

where the elasticity coefficients atj k £ are defined by 

( 1 -3-49) a ijk t = ASijS kc + 2/u6 ik S jC . 

The mapping E — > A(tr E)l + 2/uE is invertible if and only if /u(3A + 
2/u) + o (and we know that /u(3A + 2/j.) > O from above). Thus given 
there corresponds a unique Z. However, this is not always true in actual 
experiments for large deformations. This model can be expected to be 
acceptable only for small strains E. 

45 Exercises 

1.3-1. Given a matrix A € §> show that A 1 ^ 2 is uniquely defined in §>. 
If F is an invertible matrix and F = RU = VS, U = (F T F) 1 ^ 2 , 
V = (FF T ) 1/2 show that R = FU~ l ,S = V' l F are orthogonal. 
Show also that S = R, thus proving fheorem [l.3.1l 

1.3-2. If 23r is any reference configuration of an isotropic, material 
frame-indifferent material, explain why 'Lr(I) is just a multiple 
of I as shown in Theorem ll.3.61 

1.3-3. Complete the details in the proof that A,fi > O for a natural 
state. In particular, prove relations (ll.3-32t . (I1.3-34t . dl.3-37l >. 
and (ll.3-4H . 



1.4. Hyperelasticity 
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1.4 Hyperelasticity 

If the constitutive equation is taken into account, the equilibrium equa- 
tions in the reference configuration reduce to a system of three equa- 
tions, for the three components of the deformation <p, along with bound- 
ary conditions: 

(1.4-1) DIV R f R (Vcf>) + PR b R = Om& R 

(1.4-2) t R (V<f>)n R = tut on dS lR 

(1.4-3) <f> = <f> ond% 0R . 

This is equivalent to the variational equations 
(1.4-4) J f R (V<p):V6dX R = Jp R b R .6dX R + J t lR .6dA R 



LR 



for all 6 : 93 R -> R 3 , vanishing on dS oR . 

It was seen in section \L7\ that if the body forces and applied surfaces 46 
were conservative, then dl .4-4b could be written in the form 

(1.4-5) J f*(V0) : VcpdX R = B'(cf>)6 + T[(<f>)9 

for real- valued functionals B and T\ (cf. (II .2-32b and (ll.2-35t . 
If it were possible to write 

J f«(V0) : V6dX R 

'Sr 

as W (0)# for some functional W, then the problem (ll.4-4t would reduce 
to finding the stationary points of the functional W — (B + T\). 

Note that upto now, the equations which give the symmetry of X# = 
(y(f>y i T R have not been mentioned; it will be seen later (cf. Theorem 
ll.4.3t that for materials under consideration in this section, these equa- 
tions will automatically be satisfied. 

The above considerations lead to the following definition: 
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1. Description of Three - Dimensional Elasticity 



A homogeneous elastic material is said to be hyperelastic if there 
exists a differentiable function W : M 3 —> R such that 



(1-4-6) T R (F) = —(F) 

dF 

for all F € M 3 , or componentwise, 

dW 

(1-4-7) T Rij (F) = —(F). 

or ij 

A word on notation: the Frechet derivtive W(F) : M 3 — > R is a 
continuous linear operator such that for F, and F + G in M 3 , 

<W{F + G) = <W(F) + < W\F)G + o(G) 

= win + gprWGij + o(G). 

d'W 

The term - — (F)Gn will also be written as 

dFij 

dF dFij 

where the matrix — — (F) has components - — (F). 

dF dFjj 

Theorem 1.4.1. Consider a homogeneous hyperlastio material acted 
on by body and applied surface forces which are conservative. Then the 
boundary value problem with respect to <f> is formally equivalent to 

(1.4-8) f (0)0 = 

for all 6 : 23r — > K 3 , vanishing on dS () R where, for all (A : Br — > R. 3 , 
(1.4-9) m = J WiViJfiX^dXR-iBi^ + T^)). 
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Proof. Given if/ : 23r — > R 3 and 'W : M+ — > R, let 



W((A) d = I W(V4>)dX R . 



Then given tA and 8, 



Wy, + 9)- W ty) = J rWWoXXjd) - <W(Vi{,(X R )))dx R 

\—(y,MX R )) : VfHX H ) +■ "(MA' A >)|:A' A .) 



f fjt(V^r) : Vftffji + o(P||). 



Thus, at least formally, 



r f ty)6 = f 



(1.4-10) W'{ip)8 = T R (Vif/) : V8dX R 

Ik 

and the result follows. □ 

Remark 1.4.1. It must be verified in each circumstance thet W is Frechet 
differentiable and that the right hand side of ( ll.4-10t does indeed give 
the Frechet derivative. If the C 1 -uniform norm is chosen for the space of 
differentiable vector functions on S R and if the first partial derivatives 
of T Rjj are Lipschitz Continuous it can be it can be seen that is indeed 
the case. 

The functional W is called the strain energy and / is called the total 
energy. The function "W : M 3 — » BL is called the stored energy function. 

Notice that the boundary value problerm is precisely the Euler equa- 
tions associated to the total energy. 

If (f> on d$$ oR is extended to the whole of S R and I defined by 



Ioty) = I(lf/ + <f>o) 
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1. Description of Three - Dimensional Elasticity 



then one looks for <f> - <f> vanishing on <923 r such that 

f (<P - <PoW = O 

for all 6 : 23# — > R 3 vanishing on 323 r. Thus particular solutions are 
those <p which satisfy 



'inf 



(1-4-11) K</>) 



iff : SS R -» R 3 /((A). 
(A = 0o on d^ oR 



In the next chapter, it will be seen that the formulation in terms of 
the boundary value problem will be the basis for proving existence of 
solutions via the implicit function theorem while ( 11.4-111 ) will be the 
basis for the approach due to J. BALL. 

A stored energy function *W : M 3 — > R will be said to be mate- 

irw 

rial frame indifferent (resp isotropic) if Tr - is material frame 

oF 

imdifferent (resp. isotropic). 

Now necessary and sufficient conditions for a stored energy function 
to be material frame indifferent or/and isotropic will be examined. 

Theorem 1.4.2. The stored energy function 'W : M 3 — > R is material 
frame indifferent if and only if for all FeM 3 and for all QeO\ 

(1.4-12) WiQF) = <W(F). 



Equivalently, it is material frame indifferent if and only if there exists 
a function 'W : S 3 — * R such that for all FeM 3 

(1.4-13) r W(F) = <W(F T F) 

(cf. Equation dl.3-61 ). 

Proof. Since material frame indifference is equivalent to 

f(QF) = Qt{F)Q T 
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for all FeM+ and for all Qe<0\, and since 

f R (F) = det FT (F)F~ T 
it follows that material frame indifference is equivalent to 
(1.4-14) t R {QF) = Qt R {F) 

for all FeM\ and QeO\ , i.e., 

(1-4-15) —0010 - Q^(F) 

in case of hyperelastic materials. Define 
(1.4-16) 1V e (F) = WCgFJ.geO^FeM^. 

Now ifF + GeM 3 , 

W Q (F + G) = <W{QF + QG) = -^r(QF) : <2G + o(gG) 

= <2 r ^ F (!21 , ):G + (G) 
dF 

where the relation A : BC = B T C has been used (cf. Remark ["l.4.2l i. 
Thus 

&Wq T dTW 

~^(F) = Q T -^riQF). 
oF oF 

Thus material frame indifference is equivalent to 
(1.4-17) sp^Q^ ~ = °- 

□ 

Now is connected in M 3 (cf,. Exercise 1.4-2 and so the above is 51 
equivalent to) 

(1.4-18) WiQF) = <W{F) + C(Q), 
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1. Description of Three - Dimensional Elasticity 



for all FeM\, QeO\. Setting F = I,Q,Q 2 ,. . . successively, it follows 
that 

1VG2) - *W(I) + C(Q) 
<W(Q 2 ) = <W{Q) + C(G) 



Thus fo any integer p > 1, 
(1.4-19) W(<2 P ) = + PC(Q). 

Then 

W^l > P\C(Q)\ - rW)l- 

Thus if C(0 * 0, then | < W(g p )| -> +oo as p -> oo. But the set 3 + 
is compact in M+ and *W is continuous since it is differentiable. Hence 
C(Q) = and the first assertion is proved. 

To prove the second equivalence, let F = RU be the polar factoriza- 
tion of F (cf. Theorem ll.3.U . For Ce§>, set 

(1.4-20) < W(C) = <W{C 112 ). 

Then 

<W(F) = <W{RU) = <W(U) = Will 2 ) = <W(F T F) 

52 since U 2 = F T F. Conversely, if ( 11.4-131 ) is true, then for FeM+ and 
Qe0 3 + , 

<W{QF) = < W(F T Q T QF) = <W(F T F) = <W(F). 

It can be show that (cf. Exercise 1.4-4) if *W is differentiable, so is 

&W 

f W. Without loss of generality, it may be assumed that the matrix 

oC 

is symetric. For 

*W(C) = f W( C+ 1 C ), CeSl. 
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Remark 1.4.2. The following identities involving the scalar product : in 
M 3 are useful. 

(1.4-21) A : BC - ir(AC T B T ) = tr(B T ~AC T ) = B T A : C 

(1.4-22) 

A :BC = AC T :B = B: AC T = \x{BCA T ) = tr{CA T B) = CA T : B T . 

The identity (ll.4-21t was used in the proof of the above theorem. 

The following theorem says that in case of frame indifferent hy- 
perelastic materials, the symmetry of the second Piola-Kirchhoff stress 
tensor is automatically verified. 

Theorem 1.4.3. Let the material be hyperelastic and material frame 
indifferent. Then 

T d'W T 

(1.4-23) 2* = t R {F) = Y,r(F F) = 2— —(C), C = F F. 

dC 

Thus the second Piola-Kirchhoff stress tensor is automatically sym- 
metric. Conversely, if there exists a mapping <W : S 3 -> R such that 

d'W T 

(1.4-24) L*(F> = 2— (F T F) 

oC 

then the material is hyperelastic with 53 
(1.4-25) <W{F) = <W(F T F) 

and consequently is material frame indifferent. 

, - , d'W 

Proof. t R (F) = F Tr(F) = F~ l —-(F). 

oF 

Also <W(F) = <W(F T F). Now if F, F + GeM 3 + , 

<W(F + G) - WiF) = < W(F T F + F T G + G T F + G T G) - < W{F T F) 
dTW 



dC 



(F T F) : (F T G + G T F) + o(G) 



T dfW T T 

dC oC 
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1. Description of Three - Dimensional Elasticity 



by d 1.4-2 ID - dl.4-22t . But is symmetric. Thus 

dC 

&w T 

W(F + G) - <W(F) = 2F(— — (F F)) : G + o(G). 

oC 

54 Hence 



(F) - 2F"-^{F T F) 



dF dC 

or I«(F) - F" 1 —(F) - 2— — (F F). 

oF oC 

Conversely, if *W(F) = "lV(F r F) then consider the mapping F 
F r F from M+ into S 3 + . One has 

<W'(F)G = nV'(F r F)(F T G + G T F) 



or 



^rr(F) : G - ~^7(F T F) : (F G + G F) 
oF oC 



= 2F—^-{F F) : G as before. 
d<W 



dC 

Hence 

(F) = FH R (F) = T R (F) 



dF 

and the result follows. □ 

Now the effect of isotropy on a stored energy function can be simi- 
larly examined. 

Theorem 1.4.4. A stored energy function 'W : M+ — » R. isotropic if, 
and only if for every FeM.\ and for every QeG^, 

(1.4-26) <W(F) = <W{FQ). 

Proof. The argument runs along the same lines of that of Theorem ll.4.21 
and is left as an exercise (cf. Exercise 1.4-5). □ 
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Theorem 1.4.5. A stored energy function < W : M+ — > R is material 
frame indiffernt and isotropic if, and only if, there exists a function 

= Goj +°°[ ) 3 — > R 
swc/j ?/zaf = ^(; F r F ) = <P{i ff t) (1.4 - 27) 

/or every FeM+. 

Proof. By the material frame indifference, there exists a function "W : 
§> — > R such that 

<W(F) = 'W(F T F). 
By the isotropy, if QeO\, then 

- 1V(Fg) = < W(Q T F T FQ). 

55 

Thus *W : Si -> R satisfies 

<W{C) = W{Q T CQ) 

for every C £ §> and for every 2 e 0+ (since for every CeS> there 
corresponds F = C l ^ 2 e M+ with C = F T F). Now it was shown in 
the proof of the Rivlin-Ericksen Theorem (Theorem ll.3.4t that such a 
function must be a function of the principal invariants 

Cnversely, if 1V(F) - <£( W), let <2eG>+. 

Then 

l (FQ) T FQ - l Q T F T FQ ~ l F T F 
1 {QF) T QF - l F T F 

and so 

<W(F) = W{QF) = <W(FQ) 
and the thoerem is proved. □ 

The next result expresses the Piola-Kirchhoff stress tensors in terms 
of the function <p of the above theorem. 
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1. Description of Three - Dimensional Elasticity 



Theorem 1.4.6. Given a function cp : Qo, +°°[) 3 — > R and a stored 
energy function 

<W(F) = 0(n(C),i 2 (C), 13(C)), C = F r F, 

then 

(1.4-28) ^f R (F) = + |^0i/ - FF T F) + % 3 F~ T 

2 oil 012 013 

56 where 

T dcp dcp 

h = ik(F F) and — - — (i F r F ),k - 1,2,3. 
oik oi k 

Furthere 

1~ deb dcp dd> , 

(1.4-29) -X R (C) = + -%/ - C) + /^C- 1 

2 OTl 012 013 

= 1 1\ -I 12 I 

\di\ di 2 013 I 



di2 di3 j 813 

Proof Let T be the map r : M 3 — > § 3 given by F(F) = F T F. Now 

Tr(F) = —(F) where 
oF 

frW dcp 
—(F) : G = -f(i c )T'(F)G. 
OF di k 

Now i\(C) = trC and so 
(1.4-30) iv(QD = tr(D), 

13(C) = 1 [3(trC) 3 - 3(trC)tr(C 2 ) + 2tr(C 3 )] 

and so 

4(C)D = i [3(tr C) 2 (tr D) - 3(tr D) tr(C 2 ) - 6(tr C) tr(CD) + 6 tr(C 2 D)] 
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= l - [(tr C) 2 - tr(C 2 )] tr(D) + tr(C 2 D) - tr(C)tr(CD) 



= i 2 (C) tr(D) + tr(C 2 D) - tr(C) tr(CD). 



Now 



tr(C 2 D) - tr(C)tr(CD) = tr((C 2 - «i(C)C)£>) 



= \x(i 3 (C)C~ l D - i 2 D 
using the Cayley-Hamilton theorem. Thus 57 
(1.4-31) i y {C)D = i 3 {C)tr{C^D). 

Finally 

i 2 (C) = ^[(trC) 2 -tr(C 2 )] 

j 2 '(C)D - tr(C)trD - tr(CD) 
= tr((ii(C)/-C)D) 
- tr(( i2 (C)C" 1 - z 3 (QC 2 )D) 

again using the Cayley-Hamilton theorem. This may be again written as 

(1.4-32) i 2 '(C)D = i 3 (C)tr(C =1 )tr(C _1 £>) - i 3 (C) tr(<T 2 D). 

Also, F(F)G - F T G + G T F. Thus' 

dTW dd> T T 

-(F) : G = -f tv{F T G + G F) 



dF d H 



+ HC- 1 )HC~ 1 (F T G + G T F)) 

Ol2 

+ ^-h tr(C- 2 )(F T G + G T F)) 

Ol2 

+ ^-h tr(C^)(F T G + G T F)) 



Now, iv(F T G + G T F) - 2F : G 
tr(C'\F T G + G T F)) = C' T : (F T G + G T F) 
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1. Description of Three - Dimensional Elasticity 



= F~ l F~ T : (F T G + G T F) 

= 2F~ T : G(Using(QT£7i} - - HI .4-221 ) 

tr(C~ 2 (F T G + G T F)) = C~ 2T : (F T G + G T F) 
= C~ 2 : {F T G + G r F) 
= F- l F~ T F~ l F- T : (F T G + G r F) 
- IF - T F~ l F T : G 



Hence 

1#W,„ „ dO 



(F):G = ^i 3 tr(C~ l )F- T : G 



2 dF dn 



dO T i t dO T 

— i^F- T F- l F~ T : G + — i 3 F~ T : G 

Ol2 dl 3 



Now, consider 



i 3 tr (C _1 )F _1 - i 3 F- T F- l F- T 



= i 3 [tr(C~ 1 )F~ T F' 1 - F^F^'f^F'jf 
- i 3 [tr^" 1 )^" 1 - B" 2 ] F,B = FF T 



since B and C are similar. Again = ik{B). Now by the Cayley- 
Hamilton theorem, 

i 2 B~ l - i 3 B~ 2 = h I-B = ill - FF T . 

Combining all these relation dl.4-28b follows. To obtain d 1.4-291 
note that f, R {F) = F- l f R (F). Hence 

1 a dO dO T dO 1 T 

-1 R (F) = —I+ —{Hi - F F) + —i 3 F- l F- T 

2 oi\ 012 di 3 

which gives the first relation. To get the second, by the Cayley- Hamil- 
ton theorem, 



i 3 C~ l =C 2 -iiC + i 2 I 



and the result follows. 
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Remark 1.4.3. Compare the last relation in (ll.4-29t with the statement 
of Rivlin- Ericksen theorem (Theorem ll.3.4t 

Theorem 1.4.7. Consider a St venant- Kirohhoff material with 

(1.4-33) L* (£) - Atr (E) + 2/iE 

It is hyperlastic with 

(1.4-34) <W*{E) = ^(trE) 2 + fi tr(E 2 ) 

= {A+ ^\ n - 3) 2 +Kh - 3) - ^(i 2 - 3) = <Dfc) 



where 



Proof. Set 



Now 



ik = ik(C),k = 1,2,3. 
<W(C) = + 2E) = IV * (£). 



IV * (E + H) = <W * (E) + A tr E tr H + 2p tr {EH) + o(H) 
= <W*(E) + (A(trE)I + 2/jE) : H + o(H). 



Hence 



This implies that 



7? 



snV 

and hence the material is hyperelastic. The verification of the expression 
for O is left as an exercise to the reader. □ 



Remark 1.4.4. The above result gives another proof (cf. equation 60 
( 11.3-471 )) that St Venat-Kirchhoff materials are isotropic and material 
frame indifferent. 
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1. Description of Three - Dimensional Elasticity 



Remark 1.4.5. Other examples of hyperelastic materials will be seen in 
Chapter|2(Ogden's materials) 

Theorem 1.4.8. Let 23# be a natural state of a material which is isotro- 
pic and material frame indijfrent. Then if < WeC l (M.\;R) 

(1.4-35) r W*(E) = ^(tr£) 2 + ptr(E 2 ) + o(\E\ 2 ). 



Proof. Let 



Thus, 



<W*{E) = j{tr{E)f + H tr(£ 2 ) + dW * (E) 

dTW* ddW* 
—-E = A(trE)I + 2pE + -r=^(E) 
oE oE 

= 1T R (E) = AtiE) + 2pE + o(E). 



dS<W* 

-(E) = o(E). 



dE 

Since subtracting a constant in < W* does not change the stress ten- 
sors, it can be assumed, without loss of generaliry, that S'W * (o) = o. 
Hence 

i 

6<W*)(E) = J -qz-(E) : dt = o(\E\ 2 ). 
o 



61 Exercises 

1.4-1 For a non-homogeneous hyperelastic material, 

Tr(X r ,F) = —(X r ,F) 
oF 

for every Xr€Br, and for every FeM+. Extend the analyisi of this 
section to such materials. 
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1.4-2 (i) Show that M 3 is a connected subset of M 3 . 

(ii) Show by an example that M 3 is not convex. 

(iii) Identify its convex hull in M 3 . 

1.4-3 Assume that (cf. Proof of Theorem ll.4.2t 

-W(QF) = <W(F) + C{Q) 

For all FeM 3 + and Q fQ 3 . Show that C(Q) = o without using the 
continuity of "W. 

1.4-4 (i) Show that § 3 is an open subset of S 3 . 

(ii) Show that if "W is differentiable, so is "W. 

1.4-5 Prove Theorem ll.4.41 show that isotropy is equivalent to 

Tr{FQ) = t R {F)Qfov allFeM 3 

and QeO>\ , which is in turn equivalent to 

cfW d'W t d"W 

oF oF oF 

1.4-6 Check the second relation in equation dl.4-341 ). 

1.4-7 Consider an elastic matrial with 

T(B) - B (i B )I + B\{i B )B + B 2 {i B )B 2 . 

Find necessary and sufficient conditions on 62 

B a : (]o,+oo[) 3 — > R,a = 0, 1,2 

for the material to be hyperelastic. 

1.4-8 In Theorem 11.4.81 compute the terms of order 2 in Yj*r(E) and 
terms of order 3 in r lV*(£'). Explain the discrepancy in the number 
of terms obtained in each case. 
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Chapter 2 



Some Mathematical Aspects 
of Three-Dimensional 
Elasticity 

IN THIS CHAPTER, questions of existence of solutions to the bound- 
ary value of three-dimensional elasticity will be examined. In the first 
section, some general considerations about these boundary value prob- 
lems will be mentioned. The problems will be classified with respect 
to boundary conditions. As good models of elasticity must preclude 
uniqueness of solution, several examples of non-uniqueness will be pre- 
sented. 

The first tool for the study of existence of solutions is the implicit 
function theorem. As this requires a knowledge of the linearized prob- 
lem, the second section will briefly present linear elasticity and the third 
section will prove existence, albeit for a very narrow class of boundary 
conditions. The fourth section will study incremental methods, whose 
analysis follows closely related lines. 

The last two sections will present results on polyconvexity and ex- 
istence of solutions to the problem of minimizing the energy using the 
approach of J. BALL. Though several types of boundary conditions can 
be studied here, the main drawback is the lack of regularity of solutions 
and so it is not known if the solutions satisfy the equilibrium equations. 
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2. Some Mathematical Aspects of Three-Dimensional . 



2.1 General Considerations About the Boundary 
Value Problems of Three-Dimensional Elasticity 

64 

Given a respones function Tr : M 3 — > M 3 satisfying 
(2.1-1) F(f R (F)) T = f R (F)F T 

for every FeM 3 and given the denstiy p R : B R — » R and dead loads 
: Br — » I. 3 , «r : <9Si# — > R 3 , the boundary value problem arising out 
of the equilibrium equations amounts to finding a deformation CD which 
satisfies 

(2. 1-2) DIV R t R (y<f>) + p R b R = 0mB R 

(2.1-3) t R <y<f>)n R = tut on dB 1R 

(2.1-4) <p = (f> (given) on 8R oR 

The boundary condition <p = <p () on dE oR is called a boundary con- 
dition of place. The boundary condition ( I2.1-3I I on <9B 1R is called a 
boundary condition of traction (and this defnition implies it is a dead 
load). 

If dB oR - <p the problem is a /?wre traction boundary value prob- 
lem.If dB\ R the problem is a pure displacement problem. If both dB oR 
and dB\ R have strictly positive dA R - measure, then the problem is a 
mixed displacement-tranction problem. 

Recall that 

£/V0) = (V^WO) 
where ^ : M 3 —> § 3 

65 and the boundary value problem <I2. l-2b — (12. 1-41) can be rewritten in 
terms of this tensor. If the material is hyperelastic (cf. Section [L4l ). then 



(2.1-5) 



2.1. General Considerations. 



55 



for a stored energy function 'W : M 3 — > R and the problem is equivalent 
to finding the sationary points of an energy functional, 



The partial differentail equations H.l-21 are nonliner with respect to 
cp since the mapping Tr : M 3 — > M 3 is non-linear, and of the second 
order. The non-linearity occurs in the highest order terms and this is a 
source of difficulty. Another source of difficulty is that the solution <p 
must satisfy det (V0) > o. Thus for instance in (I2.1-6t ift must vary over 
M 3 which is clearly not a vector sapce; in fact it is not even a convex 
subset of M 3 (cf. Exercise 1.4-2). 

The boundary condition of traction could be replaced by the so- 
called boundary condition of pressure (which is not a dead load, but 
it is conservative). Again it is possible to have a pure pressure bound- 
ary value problem (for instance, a soap bubble or a submarine) or mixed 
displacement-pressure boundary value problems. 

These boundary condition, though being the only ones to be con- 
sidered here, are far from exhaustive. Other types of conditions are 
possible. 

In practice one can have unilateral conditions. For instance, if the 66 
body must remain in above the plane spanned by e\,e2 the boundary 
condition is > o or t/>j (Xr) > o 





Figure 2.1.1: 



It is also possible to have a mixture of displacement and stress boun- 
dary conditions. Consider Br to be a plate with a pressure p compress- 
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ing the lateral surface (fig |2.1.2l ) and given only by its horizantal aver- 
age. Then the conditions are 




Figure 2.1.2: 

67 

Apart from possibly the problem of bodies moving with constant 
velocity in a fluid, the pure traction problems are less common. The 
pure displacement problems are quite unrealistic. 

In general, several deformed states are possible for the same sys- 
tem of forces, though they may not all be physically feasible or 'sta- 
ble' . Nevertheless the mathematical model cannot recognize the feasi- 
ble or stable ones. Hence a good model will always account for non- 
uniqueness of solutions. Several examlpes of non-uniqueness will now 
be given . 

Example 2.1.1. A mixed displacement-traction problem. Consider a 
long circular cylinder fixed at either end. The body force is just its 
weight. On the lateral surface tm = o. Assume the body to be ex- 
tremely pliable, and rotate one end by an angle of 27r and reglue it in its 
original position. Then a line parallel to the axis on the lateral surface 
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will deform into a curve and thus gives another solution other than the 
natural one, which will just be a slight bending of the cylinder under 
its weight. It is theoretically possible to rotate the face by 2kn for any 
positive integer k. Thus the model must account for an infinite number 
of solutions. 




Figure 2.1.3: 

Example 2.1.2 (F. JHON). A pure displacement problem. Consider the 
body to be betwen two concentric spheres. Assume u = on both the 
inner and outer surfaces. Apart from the trivial solution, it is possible to 
have an infinite number of solutions by (theoretically) rotating the inner 
sphere about an axis by an angle of 2kn and re-glueing it to the body. 

u = 




Figure 2.1.4: 
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Example 2.1.3 (C. ERICKSEN). A pure traction problem.A rectangular 
lock is pulled normal to the upper and lower surfaces. Rotation of the 
configuration by n produces a (though urrealistic) solution where the 
body is comprssed! 




Figure 2.1.5: 



Example 2.1.4. Consider a thin circular plate subjected to the boundary 
condition d2.1-7t - d2.1-8t with p = Ap\ . If A < (i.e the plate is pulled) 
or if A > o and small, u = o is the only possible solution. If A exceeds 
a critical value, the plate can buckle upwards or downwards thus giving 
two additional solutions (cf . 12.1.61) This is a buckling phenomenon 



Example 2.1.5. Eversion problems. A cut tennis ball (borrowed from a 
very good friend) can be made to exist in two diffrent states as shown in 
fig |2.1.7l The everted state can be produced by pushing hard enough on 
the natural state. 
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Figure 2.1.6: 




Figure 2.1.7: 



It is possible to do the same thing to a tube. These are examples of 
multiple solutions to a pure traction problem. 

Returning to the various restrictions on the model and its analysis, 
the first is the taking into account of properties like isotropy, hypere- 
lasticity and the axiom of material frame indifference. These are fairly 
easy to handle. In sections [T31 and various necessary and sufficient 
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conditions of the relevant functions were studied. The main restriction 
is that the solution <p must satisfy forces 

det(V0) > o. 

In using the implicit function theorem, this requirement is ignored 
at first and then shown to be satisfied for sufficiently small forces. 
72 In the different approach of J. BALL, this is taken into account (a.e 

in Br) by imposing it on the set % of test functions over which the 
energy is minimized. This precludes the convexity of % which makes 
minimization more difficult than usual. In this approach, it will be nicely 
taken into account by imposing that < W(¥") — > +oo when det (F) — > o + . 

Even if det > o everywhere on Br, it does not ensure that (pis a 
one-one mapping, a property natural to expect in a deformation.Thus if 
a body as in fig. I2.1.8r a) in contact with the horizontal plane is pushed 
along the two 'arms', it must take a shape as in fig I2.1.8r b). But the 
mathematical model will not preclude a situation as in fig |2.1.8> c) where 
the material penetrates itself, still keeping det V0 > 




Figure 2.1.8: 

73 In the case of incompressible materials, the energy is minimized 

over a set of function \p (in a suitable function space) with 

det(VtA) = 1 a.e. 

Some of the noations used hitherto will be changed. 93« will hence- 
forth be denoted by Cl, Q. a bounded open subset of K 3 and its boundary 
(923j? by T. The portions 323 Rand<323iR will be denoted by T and T\ 
respectively. 



2.1. General Considerations. 



61 



The generic point Xr will henceforth be labelled x and dX^ and dAf> 

d 

will be changed to dx and da respectively. The derivatives - — will 

dX Ri 

be denoted by 3; and DIVr by div. The normal ur to d$$R will now be 
given by v = (v ; ), the unit outer normal to T. 

The tensors £ = (£) and Tr - (Tr.) will be denoted by (<x; 7 ) and 

R Rij 

(tij) respectively. The vectors pRbR and t\R will be denoted by / = (ft) 
and g = (gi) respectively. The symbols for cp,u,F = V0, B - FF T , C - 

T 1 

F and E - -(C — /) will remain unchanged. 

Thus, for instance the equations of equilibrium im terms of read 
in the old natation as: 



(2.1-9) DIV S 



V R ) 



+ o-RbR = o in 



(2.1-10) V<p^nR = t m ondK lR 

R 

(2.1-11) = 0oond93oR 

These, when translated in to the new notations will read as 74 

(2.1-12) -dj(o- k jd k <pi) = fi in Q. 

(2.1-13) o- kj d k (piVj = gi on Ti 

(2.1-14) 4>i = <hi on To 

Exercises 

2.1-1 . Assume that a pure traction problem has a solution <f>. Show that 
J p R b R dX R + J t lR dA R = 

and J" (pKpRbRdXR + J" (pAt[RdAR = 0. 
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75 



2.1-2 .Consider a hyperelastic incompressible material. In the con- 
strained minimization problem 



where 



inf /((A) 
<U = {ift I det(VtA) = la.e.}, 



show by a formal computation that the Lagrange multiplier is the 
pressure. 

2.2 The Linearized System of Elasticity 



Consider the boundary value problem d2.1-12l )- (l2.1-14t . In terms of the 
displacement u it can be rewritten as 



(2.2-1) 
(2.2-2) 
(2.2-3) 



-djiq-ij + cr k jd k Ui) = ft in Q, 
(o-ij + cr k jd k Ui)vj = gi on T u 
u t = u oi on T , 



with the constitutive equation 

(2.2-4) o-ij = o-*j(E(u)) = AE kk {u)5ij + 2fiEij(u) + o{E) 

where 

(2.2-5) E(u) = - (Vu T + Vu + Vu T Vu) 

If u were defined in a suitable function space, whose functions van- 
ish on To, then symbolically one can write 



(2.2-6) 



A{u) = 



f 



The linearised system of elasticity will then be formally defined as 
(assuming A is differentiable) 



A'(0)w 
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This can be derived as follows. The linearized strain tensor is 
(2.2-7) e{u) = - (Vu T + Vw) . 

Then o~ will in turn be linearized as 
(2.2-8) (Tij = Ae kk (u)6ij + 2fieij(u). 

Substituting this in ( 12.2-1 t - (!2.2-3t and keeping only the first order 
terms, the linearized syatem elasticity turns out to be 

(2.2-9) -djo-ij =fi in Q 

(2.2-10) o-ijVj = gj onTi 

(2.2-11) Uj - u j on To 

where cr is given by ( I2.2-8I ). Note that such a syatem cannot be a model 
for elasticity (of Exercise 2.2-1) but only approximation of a model. 

Remark 2.2.1. If the equations were written in terms of tij and then 
linearized, the same linearized system of elasticity would have been ob- 
tained. This is because Tr = (I + Vm) £ and on linearizing this realtion 

R 

only the part coming from / 2 = Z Wli l be retained. 

R R 

Before the existence and regularity of solution to the linearized sys- 
tem of elasticity can be studied the following notations for the Sobolev 
spaces will be needed. 

Let m > be an integer and 1 < p < +oo. Then 

(2.2-12) W" h '\Q) = {veL p (Q.) \ d a veL p {Q)iox all | a \< m) 

where a is a multi-index and d a v is the corresponding partial derivative 
( in the sense of distribution). This space is a Banach space with the 
norm 



(2.2-13) || v \\ m , p . n -- 



f ^ | d a v f dx 
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(with the standrad modification if p = +00). The semi-norm | . L^n is 
defined by 

( \ l 'P 
J] I d°v \ p dx 



(2.2-14) || v \\ m , n = 



If m = 0, W°-P(fl) = Z/(£2)and | . | 0jA q is the usual LP {ay norm. 
If i^(f2) is the space of C°° functions with compact support in Q., its 
closure in W m,p (Q) will be denoted by W™' p (Cl). 

If p = 2, it is customary to write H' n {Cl)m&H™(Q.) instead of W m , 
2(Q) and V^™ ,2 (Q) respectively. The norms and semi-norms in this case 
will be written as || • || m> £> respectively | • | Wj n is the L 2 (Q.)- norm. 

By Poincares inequality, | • \ m , p £i is a norm on W™' P {Q.) and is equiv- 
alent to || . || m , p n,forl < p < 00. 

In case of vector valued or tensor valued functions, the symbols 
W m ' p (Q.), L p (0) will be used to denote that each component is in W m,p 
(Q.) or L P (Q.) respectively. However the symbols for the norms and semi- 
norms will not be altered. 

The following result is fundamental. 

Theorem 2.2.1 (Korn's Inequality). Let Y be smooth enough. Then 

(2.2-15) {v = (v,-)eC 2 (0) | e, 7 (v)eL 2 (0), 1 < i, j < 3} = M l (Q) 

Consequently, there exists cons tans C\ > 0andC2 > such that 
(2.2-16) 

Ci || v ||!, n < (| v | 2 n + I 6(v) | 2 n ) 1/2 < C 2 || v \\ UC iforallveH\a). 

Proof. See DUVAUT and LIONS [1972] or NITSCHE [1981]. The 
main difficulty is in proving d2.2-15t . Since the second inquality of 
( 12.2-161 ) is obvious, the first follows from ( 12.2-151 ) and the closed graph 
theorem. □ 

A consequence of the above result is 

Theorem 2.2.2. Let V be defined by 

(2.2-17) V = (veH^O) | v = on T } 
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where the da-measure of To is strictly positive. Then the semi-norm 
I € ( v ) \o,n is a norm on V equivalent to the norm \\ . \\i t a. 

Proof. Cf. Exercise 2.2-2. □ 

Assume now that u = on To. Let V be as in d2.2-17t . Multiplying 
( 12.2-91 by a function veV, integration by parts using Green's formula, 
and using (12.2- 10b . (12.2-llt and the symmetry of o~, the following vari- 
ational formulation of the problem d2.2-9t - (I2.2-1H can be obtained. 

Find ueY such that, for all veY, 

(2.2-18) a{u,v) = L{v) 

where 

(2.2-19) a(u,v) = J (Ae kk (u)e a (v) + 2^e,7(«)ey(v)) 



and 

(2.2-20) L(v) - fvidX + J gm da. 

n r. 

By a simple application of Schwarz's inequality, it follows that a(., .) 
is a continuous bilinear form on H (H) and L is a continuous functional 
on H'(£2) (and hence on V as well). 

The following existence result holds. 

Theorem 2.2.3. Consider the variatiational formulation of the linea- 
rized syatem of elasticity, (I2.2-18I I. or, equivalently, the problem: Find 
ueY such that 

(2.2-21) J(u) = infJ(v) 

veV 

where 

(2.2-22) J{v) = ^a(v, v) - L(v) 

ifA>0 and p > then the problem has a unique solution. 
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Proof. Observe that by Theorem 12. 2. 2l for all veV, 



(2.2-23) 



a(v,v) >2n\ e(v) \l <a > a \\ v \\{ 



1,Q 



where a > o (using A > o and p > o). Thus J : V — > R is a convex, and 
continuous functional. Hence it is weakly lower semi continuous. Let 
{iik) be a minimizing sequence in V, i.e. 



Since J is coercive (i.e., J{v) — > ooas || v ||— > ±oo) it follows that 
80 is a bounded sequence and hence has a weakly convergent subsequence. 
Let u be the weakly limit of the subsequence (again indexed by k for 
convenience). Then by the weak lower semi-continuity of J. 

inf J(v) < J{u) < lim inf J(uk) - inf 7(v). 

veV k veV 

Thus J attains its mimimum at u. It is easy to see that equations 
(12.2-181) are simply equivalent to the equation J'(u) = 0. Hence the 
equivalencve of the two problems since J is convex and the existence of 
a solution. 

If u\ and ii2 are solutions in V then a{u\ - U2,v) = for all veV. 
Setting v = ui - ii2 and using d2.2-23t . it follows that u\ = 112, thus 
proving the uniqueness of the solution. □ 

Remark 2.2.2. The existence of a unique solution to ( 12.2-1 8t also fol- 
lows directly from (12.2-231) by applying the Lax-Milgram Lemma. 

Finally, let us state the result on the regularity of solutions to the 
above problem. 

Theorem 2.2.4. Suppose T is smooth enough and for some p > 2, 
feC p (Q). Assume Ti - (p. Then the solution ueY - H[ ( (Q) of the corre- 
sponding linearized system of elasticity is in the space Y P (Q), where 



J{uu) — > inf J{v) < +00. 



(2.2-24) 



V p (0) - {veW 



2 ' P (H) I v - on T}. 
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Proof. The case p = 2 has been proved by NECAS [1967]. If p > 2, 
the argument goes as follows. Let A'(0) : V p (f2) — > C p (fl) repre- 
sent the operator of the linearized system of elasticity. Then, if in- 

Hef 

dex (A'(0)) = dim(kerA'(0)) - dim(CoJfcerA'(0)), it was proved by 
GEYMONAT [1965] that for all p, 1 < p < oo, the index was inde- 81 
pendent of p. Now, by the result of Necas above if p - 2,A'(0) is 
onto and so dim(Coker(A' (o)) = 0. By uniqueness of the solution, 
dim(ker(A' (0)) = for all p > 2. Hence the index is zero for all p 
and so dim(Coker(A' (o))) = for all p > 2, i.e., A'(0) is onto for all 
p > 2, which proves the theorem. □ 

Remark 2.2.3. The above result does not follow from those of AG- 
MON, DOUGLIS and NIRENBERG [1964]. Their results state that if 
feh p {Q) implise ueW 2 >P(Q.) then feW m 'P(Q.) implies ueW z+m *(Ci) for 
our problem. The 'starting' regularity result (m = 0) needs be know a 
priori and Theorem l2.2.4| proves that in case of the linearized syatem of 
elasticity. 

Caution! The W 2 ' p regularity does not hold for the mixed displace- 
ment-traction linearized syatem of elasticity. 



Exercises 

1 T 

2.2-1 Show that if e(u) = + ^ u ^' tnen a constitutive equation 

of the form cr = cr*(e), with o~* a linear function of e, does not 
satisfy the axion of material frame indenfference. 

2.2-2 Prove Theorem 2.2 - 2. (Show first that |e(v)|on) is a norm on V. 
Prove the equivalence of norms by contradiction: assume there 
exists a sequence v*W with ||v*||i,jj - 1 and |e(v A ')|o,n —> 0. 

2.2-3 Consider the linearized system of elasticity in variational form 
with Tq = (p. Show that there exists a solution to the problem 
provided 

f f dx = f dx, 
n r 
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82 in the quotient space H Y (Q.)W, where 

W = {veH^Q) | e(v) - 0}. 

Show also that 

W = {velH^miv - a + b A ox}. 

2.2-4 Extend the regularity result to the case 1 < p < 2. 

2.2-5 Show that if p > 0, there exists a Aq < such that if Ao < A < 0, 
the linear form a{., .) given by ( 12.2-191 is coercive. 

2.3 Existence Theorems via Implicit Function The- 
orem 

In this section, existence solutions to the pure displacement problem will 
be proved using the implict function theorem. 

For simplicity, consider first a St Venabt-Kirchhoff marerial Recall 
that the constitutive equation in this case can be written as 

(2.3-1) o-jj = Oijkt E k( {u) = AE kk {u)5ij + 2pEij{u), 

with A > and p > 0. Also 

(2.3-2) Eij{u) = eij{u) + -diU m djU m , 

where 

(2.3-3) £ij(u) = ^(djUj + djUi) 

Then the boundary value problem (12.2-lb — (12.2-31) becomes (T\ = <p), 

(2.3-4) - dj(a ijpq e pq + -a ajpq d p u m d q u m + a kjpq d p u q d kUi 

+ -jakjpqdpUmdqUmdkUi) = fj in Q,, 
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(2.3-5) u = on T. 

This can be written as 
(2.3-6) A{u) = f in Q, 

(2.3-7) u = on T, 

where A(u) - (Ai(u)) and 

(2.3-8) Ai(u) - -dj(a ijpq e pq (u) + ja ijpq d p u m d q u m + a kjpq d p u q d k Ui 

The following existence result holds. 

Theorem 2.3.1. Assume that T is smooth enough. Then for each p > 3 
there exist a neighbourhood & p ofO in IfQ. and a neighbourhood < U P 
ofO in 

V p (0) = {veW 2 ' p (Q)|v = on T} 

such that for every fe^ p the boundary value problem ( 12.3-61 ) - d2.3-7l l 
has one, and only one, solution in ( U P . 

Proof. Since Q. c K 3 , if p > 3 the inclusion 

W l ' p (Q) -> C°(Q) 

is continuous. Further W 1>P (Q) is an algebra (cf. ADAMS [1975]). Thus 
if f,geW l > p (Q.),fgeW Up (£l) and 

(2-3-9) \\fg\h, P ,n < C\\f\h, p n\\g\\u P .£i 

Hence A : V*(Q) c W 2 ' P (Q.) -> C P (Q) is well-defined and is in- 
finitely Frechet differentiable. (In fact D 4 A = 0). Since A(0) = 0, the 
conclusions of the theorem will stand proved if it is shown that 

A'(0)el som(V p (Q),I/(Q)) 
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by virtue of the implicit function theorem. But the problem 

(2.3-10) A'(0)k = /, ueY p (Q.) 

is none other that the linearized syatem of elasticity: 

(2.3-1 1) dja ijpq e pq (u) = fi in CI 

(2.3-12) u = on T 

A'(0) is continuous. It is one-one since the solution of the above system 
unique for p > 2. Also by the regularity theorem (cf. Theorem l2.2.4l ) it 
is onto as well. Hence by the closed graph theorem A'(0) is an isomor- 
phism form from V P (£T) in to L p (f2) and the theorem is proved. □ 

Remark 2.3.1. This proof breaks down in the case of the mixed-dis- 
placement traction problem because of the lack of W 2 ' p (f2) regularity of 
the associated linearized syatem. 

85 Remark 2.3.2. One could think of solving the problem by defining A on 
W l ' q (Q.) taking values in W 1 ' q (Q.) , thus avoiding the need of the regu- 
larity theorem which is not valid for other boundary conditions. Unfor- 
tunately, it has been proved by VALENT [1979] that on such spaces A 
is not Frechet differentiable. 

Remark 2.3.3. If were replaced by smooth functions a/^Ot) , the 
result is still true, thus extending the result to non-homogeneous mate- 
rials. 

In case of St Venant-Kirchhoff materials, it turned out that if u € 
W 2 ' p (£l), then E(u) e W l ' p (Q.) and since cr* was linear in E, cr*{E(u)) e 
W l ' p (Q.). For more general constitutive equations given cr* : S 3 — » S 3 , 
it must first be proved that if E e W 1 ' P (Q.), then 

rlef 

o-*(E)(x) = o-*(E(x)),xeQ, 

is indeed in W 1,p (0). The following result answers this question. It is 
due to VALENT [1979]. 
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Theorem 2.3.2. Let p > 3. Given a tensor field E € W l ' q (Q.) and a 
mapping cr* e C'(M \ M 3 ), the matrix valued function 

cr*{E):XeQ.^ cr*{E{x)) 

is also in W 1,9 (Q) and 

(2.3-13) d q (cr*AE(x)) = -^i(E(x))d q E k( (x). 

If cr* is of class C m+l ,m > 0, them the mapping cr* : W 1,P (Q) — > 
W 1,p (fl) defined above is of class C m and it is bounded in the sense 

(2.3-14) sup < r lltfV*(*)ll<+°° 

\\E\\i, p ,n 

for every r > 0. 86 

Proof Step®. Letcr* € C'(M\M\). Let E e W 1 ^). Then the 
components of E are all continuous and so 

cr*(E(x)) e C°(n;M 3 ) ^ L p (£2). 

Assume now that (I2.3-I3t has been proved. Then as 

OCT*: 

■(E{x)) € C°{Cl)andd q E M (x) e Z/(Q). 



it follows that d q cr*{E) e ][/(Q). Hence cr*(£) € W 1 ^(Q). 

Now d2.3-13l > will be proved. It must be shown that for any <p e 

(2.3-15) J cr* j {E{x))d q <p{x)dx = - J -^-{E(x))d q E u {x)^x)dx 

Let <p € be fixed. If £ € C^M 3 ) then I2.3-I5i follows 

by a direct application of Green's formula for smooth functions. But 
C^OjM 3 ) is dense in W l -P(Q). Thus given E e W^Q), let E„ € 
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C\D.\ M 3 ) such that E n -> E in W 1,p (fi). The relation <I2.3-15I > is valid 
for each £"„. 

Now E n — > £ in C (O : M 3 ) as well, i.e. uniformly. Thus 
are all uniformly bounded in Q and so by Lebesuge's Dominated cover- 
gence theorem. 

(2.3-16) ^o-*j{E n {x))d q <p(x)dx^> ^ cr*j(E(x))d q <p{x)dx. 



n 



Now. 

dcr* Oct* 

--±{E n {x))<p{x) —±(E(x)Mx) 

3E E c dE k 

87 uniformly and hence in L p '(Q.),p' the conjugate exponent of p. Since 
dq(E n )ke — > hi L P (Q), it follows that 
(2.3-17) 

f 5(7 t/ f 5o 1/ 

-r^{E n {x))d q {E n ) kc {x)(plx) ——(E(x))d q E kf (x)dx 
Jn dE kl J n 3E kc 

and thus (12.3-151) is established for E e W l ' p (Q.). 
Step(ii). It will be now shown that cr : W^Q) -> W 1 ' P (H) is contin- 
uous and bounded. Let E n ^> E in W 1 ' p (0). Then as before E n — > £ 



in C°(0;W :i ) as well. Hence o-* r {E n ) —> cr*..{E) uniformly and also in 
L p (Q). Similarly 

dcr*. do-* u 

— He") - —±(E) 

8E ke 8E kf 
uniformly and d q {E n ) k f — > d^-E^ in L P (Q). 'Thus by ( I2.3-13I I. 

d q (o-*j(E n )) d q (.o-*j(E)) in Z/(H), 

thereby proving that cr*(E n ) — > cr*(E) in W 1,p (Q). Thus the mapping is 
continuous. 

If Plli, P ,n < r then |£| 0>00 ,n < C(r). It then follows that cr\E) 
is bounded uniformly and hence in L P (Q) by a constant (which is a 

function r). Again (£) is bounded uniformly by a constant and 

8E kt 
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d q Eki is bounded in L P {Q) by a constant which depends only on r. These 
obsersvations lead us to the relation 

(2.3-18) sup < r ll<r * (£)lll ^ <+0 ° 

\\E\\ Up ,n 

for every r > 0. 

Step(iii). Let cr* € C 2 (M 3 ;M 3 ). It will be shown that a* : W Up (Q) is 88 
of class C 1 and that 

(2.3-19) D<rUE)G = ——^-(E)Gke- 

3 dE ke 

for any G e W L P{Q.). 

d(r*j 

By step (i), as is in C , it follows that for 

oEu 

do-*. 

E e W hp (fi), —HE) € W lp (Q). 
oE u 

Since W l ' p (Q.) is an algebra 
do-*. 

^ L {E)G M € W 1 '^) 
oE ke 

and hence the right hand side of ( 12.3-191 ) defines a continuous linear 
operator on W l ' p (£l). To show that it does indeed define the Frechet 
derivative , consider for x e Q fixed, 

(cr * tJ (E + G)- crUE) - —l{E)G ke ){x) 
1 oEu 

l 



— -{E + tG){x) -(£)(*) I dt. 

dE u 8E U " 



o 

For (£,G)eM 3 xM 3 ,let 
l 



^(E,G)^ 



/(da* a do~*ij \ 

(^7 <£ + ,G, -^7 (£) )^ 
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The mapping e¥: M 3 xM 3 ^R defined in this fashion is of class 
C 1 , since o~* is now assumed to be of class C 2 . Thus by the result of 
steps (i) and (ii), the associated mapping 

e*f: (E,G) e (W hp (Q) x W^CQ)) ->e*f (E,G) e W 1>p (n) 

is well- defined and continuouis, so that in particular , for a fixed E € 

€%(E,G) -»e /y (£,0)=0 
89 in W 1,p (n) as G ^ in W^Q). Since 

do-*. 

cr*.(£ + G) - cr*.(£) - Q^( E )Gkt = G u e» (£, G), 
it has thus been proved that 

3cr* J; ; 

Oct .(£) - _ii(£)G«. 
The continutiy of Dcr*. follows form that of the partial derivatives 

do-*. 

g^-, which is again a consequence of step (ii). 

Step (iv). To show the boundedness of Dcr * (E). Now, 

dcr* 

\\D(tUE)\\= sup \\—L(E)G ke h, P v 

\\Gh, p , a <\ oE kt 

which is readily seen to be bounded by a constant depending on r where 
l|£lli,/),fi ^ r. Thus it follows that 

(2.3-20) sup \\Do- * (E)\\ < +oo 

\\E\h,„,n<r 

for every r > 0. 

The assertions for cr* € C m+l (M 3 ; M 3 ) follow by iterating the above 
arguments. □ 
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Let u e VP(Q). The A : V*(Q) -> is defined by 

(2.3-21) A(m) = - 4- div ((/ + Vm)ct * (£(«))). 

That this indeed maps W 2 ' p (£l) into L P (Q) is a consequence of The- 
orem l2~3~2l If u e W 2 - p (Q), then Vw e W Up (Q)E(u) € W^Q) 
(since this space is an algebra). Now by the above mentioned theorem, 
<t * (E(u)) € W 1,P (Q) and as it is an algebra, (/ + Vw)cr * (£(«)) is in 90 
W 1 ' P (Q.) and its divergence is in L P (H). It is also as regular as the map 
cr* : W l ' p (Q.) — > W l,p (Q.) as the other mapping found in the map A are 
linear or bilinear. 

The boundary value problem for the pure displacement problem re- 
duces to: given / € L P (Q), find u e V P (Q) such that 

(2.3-22) A(u) - /. 

Theorem 2.3.3. Let T be smooth enough, (i) Let p > 3 and cr* e 

C 2 (M 3 ,M 3 ). ThenAmapW 2 < p (Q)intoU(Q)andisofclassC l . If 
in addition 

(2.3-23) cr*(E) = A tr(£)7 + 2ju£ + o(£) 

witfi A> and p > 0, then A\o) e Isom(Y p (Q.),h p (Q.)). 

(ii) If a* e C 3 (M 3 ,M 3 ) and ifA'(0) € from (V^ (Q), 1/(0)), then 
there exists p P > such that for all < p < p p and for all 

ve< = {veV p (Q)|||v|| 2 ,p, n < p } 
A'(v) € /ram (V P (Q), L P (Q)). Fwrffer 
(2.3-24) y P = f suplKA'Cv^l^+oo. 

A/50, ?/ze map v — » (A'(v)) -1 itf Lipschits continuous on Bp i.e., 

.Pdef ||(A-(v))- 1 -(A^(>v))- 1 || _ 
(2.3-25) L£ - sup < +00. 
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Proof. (i) That A maps W 2 - p (Q.) in L p (0) and is of class C 1 follows 
from obsevations made above. A simple computation shows that 

/Vv|Vv r \ 

A;(o)v = -dj(Do-*j(o)[ J + o-*.(o)<9*v,-). 

If o~* is of the form (I2.3-23I) . this reduces to 

(2.3-26) a;(o)v - € tt (v)<5, 7 + 2// e, 7 (v)) 

since cr*(o) = o and 

Dcr*j(o)G = AG kk 5ij + l^Gij 

But d2.3-26t is just the linearized system of elasticty (cf. Section 
12.2b and by Theorem 12.2.31 and 12.2.41 is an isomorphism as was 
shown in Theorem l2.3.1l 

(ii) Letcr* e C 3 (M 3 ,M 3 ). Then A € C 2 (W 2 -P(Q); U{Q)). Since all 
mappings occurring in A have bounded second derivatives, 

Hpf 

M p (p) = sup ||A"(v)|| < +oo. 

Il v ll2,p,£2<p 

Note that M p {p) is a non-decreasing funding of p. Now, 

(2.3-27) sup \\A'(v)-A'(o)\ <pM p (p). 

IMb,p,n<p 

If v € B p p , then 

(2.3-28) A\v) = A'(o)[I + (A'(o)y\A'(v) - A'(o))] 

But 

(2.3-29) \\(A'(o))~ l (A'(v) - A'(o))\\ < y p oP M p (p) 

where 

(2.3-30) yl = \\{A\o)T\ 
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lfp<p p where p p Q is such that 
(2.3-31) P?M p (p p ) < (yf)" 1 

then it follows from <!2.3-28b and ( 12.3-291 ) that A'(v) is an isomorphism. 
Further 



ip»)-'ii< 

1 - y^pM F {p) 



where 



(2.3-31) v" 



P def Jp 

' ~ \-y p pM p (p) 



Finally if v, w e B p , then 

(A'(v))" 1 - (A'(w))" 1 = (A'(v)y l (A'(w)-A'(v))(A'(w))- 1 
and (12.3-251) follows with 
(2.3-33) L p ^ M p (p)(y p p f. 



The latter part of the above theorem will be needed in the study of 
incremental methods (cf . Section 12.4ft . The former part leads directly 
to the following existence theorem. 

Theorem 2.3.4. Let T be smooth enough and a* e C 2 (M\M\). Let 
cr*(E) be as in (2.3-23) with A > and p > 0. Then for any p > 3, 
there exists neighbourhoods J(f p of in h p (Q) and ^ p ofO in V P (Q) 
such that for each f € <%f p there exists one and only solution u e ^ p to 
equation (!2.3-22t . 

Proof. By the previous theorem , A'{o) is an isomorphism and the result 
follows, as in Theorem l2.3.1l from the implicit function theorem. □ 



Remark 2.3.4. Theorem l2.3.1l is contained in Theorem l2.3.4l 
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The following result compares the solution as guaranteed by the 
above theorem and the solution of the linearized problem. It will thus 
be seen that for 'small' forces the linearized system is indeed a good 
approximation of the original model. 

Theorem 2.3.5. Let the assumptions of the previous theorem hold with 
cr* £ C 3 M 3 ;M 3 . For f € & p c h p (Q),letu(f) € W c V p (0) denote 
the solution to the problem (I2.3-22I) . Let ua n {f) e V P (H) denote the 
solution of the equation 

(2.3-34) A\o)uUf) = f- 

Then 

(2.3-35) MJ) - u m (f)h, P ,ci = 0(\f\l pS1 ). 

Proof. By the implicit function theorem, it follows that u is also differ- 
entialbe as a functions of /. Thus 

A'(u(f))u'(f)= I mL p (0), 

In particular, taking / = 0, it follows that 

(2.3-36) vt \o) = {A! \o))~ l . 

Now 

u(f) = u(o) + u'(o)f + o(\f\ 2 opa ) 

as A is of class C 2 . But u(o) = o. The result now follows from <!2.3-34t 
and ( 12.3-361 ). □ 

94 A major open problem in elasticity is to prove the existence of a 

solution 'close to zero' when / is 'small' , for the mixed displacement- 
traction problem. One could then compare the solutions of A(u) = f 
and A'(o)u = f. 

It was remarked in the beginning of this chapter (cf. Section 12. Il l 
that even if we solved the problem, the solution must further satisfy the 
condition det(V0) > o, and in additionm be one-one. Hitherto these cri- 
teria have been ignored. The following result assures that if / is 'small 
enough' then these conditions are met. 
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Theorem 2.3.6. Let the assumptions of theorem 12. 3.4\ hold. Further 
let r be connected. Then if \f\ . P ,n is sufficiently small, the mapping 
(j) = Id + u satisfies det (V0) > o and is one-one. 

Proof. If \f\ ,p,a is small, then \\u(f)\\ 2tP ci is small. Since W 2,p {D.) «-» 
C l (Ci), it follows that ||«||i >0 o,n is small. Hence, as the determinant is a 
continuous function of components of a matrix, it follows that 

det(V0)(jt) = det(7 + Vm)(jc) > 0, for all x € Q 

Since (f> e W 2 ' P (Q.), it can be extended to a function <f> e W 2 ' P (G) 
C\ff), where G D Q (cf. NECAS [1967]). Now (p = Id, which is 
one-one, on T. If follows from a result of DE LA VALLEE POUSSIN 
or MEISTERS and OLECH [1963]((cf. Remark lT3~3l below) that <f> is 
one-one Cl. □ 

Remark 2.3.5. The result of MEISTERS and OLECH states that if <p e 
C\G; !*), c R" an open subset, if K is a compact subset of G with 
dK connected and, finally if (p is such that det(V0) > on K and (p is 95 
one-one on dK, then <p is one-one on dK. This result can be strengthened 

o 

by allowing det V(p(x) = o on a finite subset of & and an infinite proper 
subset of dK. 

Remark 2.3.6. It is not quite necessary to resort to the use of the fairly 
deep result of MEISTERS and OLECH. If \f\ , P ,n is small ||w||i )00 ,o 
small, so without loss of generality it can be assumed that ||Vm(x)|| < 1 
for some matrix norm induced by a vector norm, for all x e Q. Now if 
4>eC°(Q.)r\C\a) and if Q. is convex, then 

\\(f>(xi) - <p(x 2 ) - (xi - x 2 )\\ = \\u(xt) - u(x\) - u(x 2 )\\ 

< sup ||Vm(x)||||xi - x 2 || 

xe]x u x 2 [ 

< 11*1 - *2||. 

Thus if x\ ± x 2 , then necessarily (f>(xi) + <p(x 2 )- 
If u is 'small' then the strain tensor 



E = -{Vu T + Vu + Vu T Vu) 
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is also 'small'. An open prolem is to study under what sufficient condi- 
tions the converse is true. 



Exercises 



2.3-1 Prove the analoge of Theorem 12. 3 . 41 with the condition u = on 
T repalced by u = uq on T. 

2.3-2 Prove the analouge of Theorem 12.3.41 using the function spaces 
C m a instead of the spaces W m -"(fl). 

2.3-3 Prove the result of MEISTERS and OLECH using the topological 
degre. Show also that in Theorem 12. 3. 61 d>(Q. = Q.). 

96 2.3-4 Give a counter example to the result of MEISERS and OLECH 
when dk is not connected. 



2.3-5 If ii > 0, show that there exists Aq < o such that if Aq < A < 0, the 
existence result of Theorem 12. 3.41 still holds. 

2.3-5 (LE DRET (1982)) Examine the existence of a solution to the 
pure displacement problem in the incompressible case, i.e., det(/+ 
Vw) - 1. 



2.4 Convergence of Semi-Discrete Incremental 
Methods 

Consider again the pure displacement problem: 

(2.4-1) - div ((/|V M )cr* (£(«)) = fin Q 

(2.4-2) u = o on T. 

It was shown in Section l231 that for small forces /, the problem had 
at least one solution which was also small. Considering that approach 
via the implicit function theorem as a 'direct' approach to the existence 
theory, by contrast the incremental methods provide a ' constructive' 
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approach to the same. Because of its constructive nature, it could be of 
use for numerical approximation of the solution. 
The basic idea is the following. Let 

< A < A 1 < . . . < A n < A n+l < . . . < A N = 1 

be a partition of the interval [0,1]. Let / be a given sufficiently small 
force in ^ p (cf. Theorem 12 .3 .41 ) . Let U" be the solution of 

(2.4-3) A(U n ) = A n f, 

U" belonging to *U? since A n f e if fe&P and JW is a ball. Note 97 
U° - o. Let u" be an approximation of U". The idea is to construct w" +1 
knowing u n , via a simpler problem, namely a linear problem. Now 

A(U n+l ) - A{U n ) = (A n+1 - A n )f. 

If A(U n+l ) is expanded about A(U n ), 

A(U n+1 ) = A(U") + A'(U n )(U n+1 - U n ) + o(\U n+1 - U"\). 

This inspires the equations defining the approximations u n . Thus 
one tries to solve the sequence of problems 

(2.4-4) A{u n ){u n+l - u n ) = {A n+1 A n )f, < n < N - 1, 

(2.4-5) u° = o. 

Of course, it is necessary that at each stage A'(u n ) be an isomor- 
phism from Y p (Q) onto LP (SI). 

The following simple, yet crucial, observation is basic to the analy- 
sis of the above method. The equations d2.4-4t - d2.4-5t can be rewritten 
as 

(2.4-6) , = (A'(u n )) f 

(2.4-7) u° = o 

which is none other than Euler's method for approximating the ordinary 
differential equation 

(2.4-8) u'(X) = (A'(u(AW l f, u(o) - o. 
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Theorem 2.4.1. Let cr* be of class C 3 (M 3 ,M 3 ). Let p > 3 and o < cr < 
Po (cf. Theorem \2.3.3\ . Let feh p (Q.) be such that 

(2-4-9) \f\ , p , a < piyPy 1 

Then the ordinary differential equation (l2.4-8t for o < A < 1 has a 
unique solution u{A) in the ball B p p . Besides 

(2.4-10) u(A) = u(Af). 

Proof. The existence of a unique solution of the ordinary differential 
equation is classical. It is converted into an integral equation and using 
the estimates of Theorem 12.3.31 regarding the uniform houndedness of 
(A'(v)) -1 and the Lipschitz continuity of the map v — > (A'(v))" 1 on Bp, 
the result follows by a use of the contraction mapping theorem. 
Now, 

u\A) = (A'(u(AW l f 

or 

{A\u{A)))u\A) = f 

or, again 

4t(A(u(A)) - Af) = 0. 
dA 

Thus 

(A(u(A)) = Af + C, 
and as w(0) = 0, C — 0. This proves the theorem. □ 

Remark 2.4.1. The equation A(u) = f has been imbedded in a one pa- 
rameter family of problems (A(u(A)) = Af, where w(l) = u. Knowing 
a solution for one value of A, here, A - 0, u(0) = 0, one tries to go 
continuously to A = 1 . This is the basis of the so - called continuation 
methods, (cf. RHEINBOLDT [1974]). 

Remark 2.4.2. The condition d2.4-9t makes precise the neighbourhood 
^ p of o in Jz? p (Q) for which a solution 'close to zero' was guaranteed 
by Theorem l2.3.4l 
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The following theorem proves the convergence of the incremental 
method desribed above when the 'mash size' max (A - A") ap- 

0<n<iV-l 

preaches zero. The assumption that the applied forces should be small 
enough is also corroborated by numerical evidence: Otherwise it is often 
observed that the approximate solutions 'blow up' for a certain critical 
value of the parameter, corresponding for example to a phenomenon of 
'buckling'. 

Theorem 2.4.2. Let o-*eC 3 (M 3 , U 3 ). Let p > 3 and < p < p p Q . If 

\f\o.p,n<p(7 P p r\ 

then given any partition 

o = A < A 1 < ... < A N = \ 
of[0, 1], Euler's method; Find u n ,0 < n < N, u n eB p , such that 
(2.4-11) A'(u n )(u n+1 - u n ) = (A n+l - A n )f 

with u° = o, is well defined and 

(2.4-12) max \\u n - u(A n f)\\ 2pQ . < C max (A n+1 A n ) 

0<n<N ' 0<n<N-\ 

where 

C = C(p,P,\f\o,p,n)>0. 

Proof. The proof is classical and again relies on the uniform bounded- 
ness and the Lipschitz continuity of the map V — > (A'(V)) -1 on Bp. 
If AA = max (A n+l - A n ), then 

0<n<N-\ 

(2.4-13) \\u N - u(f)h, P ,n = 0(AA). 

If cr N and cr were defined by 
(2.4-14) o- N = cr*{E{u N )) 
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and 



(2.4-15) 



cr = (T\E{u(f))) 



then it is easy to see that 



(2.4-16) 



To conclude, two open problems will now be stated. 

The first is to analyse the fully discrete incremental method by add- 
ing the effect of finite element methods. 

Secondly, one can construct formally a semi - discrete (or fully dis- 
crete) incremental method for the mixed displacement-traction problem. 
If it can be shown that the approximants u" exist uniquely and that they 
converge in some sense as AA — > 0, this could provide a valuable exis- 
tence theorem for this class of problems. □ 

ioi Exercises 

2.4-1 (DESTUYNDER AND GALBE (1978)). For a St Venant - Kirch- 
hoff meterial show that the map A — > u{A) (cf. equation (12.4- lOt ) 
is analytic in a neighbourhood of 0. 

2.4-2 Apply Newton's method to the equation A{u) = f, ueY p (Q.) and 
study its convergence to a solution of the equation 

2.5 An Existence Theorem for Minimizing Functio- 
nal and Outline of its Application to Nonlinear 
Elasticity 

In this section, an existence theorem for minimizing a functional will 
be proved. The functional considered will resemble the total energy 
functional of elasticity described in Section [T!4l Unf ornately, however, 
the energy functionals of elasticity will not satisfy all the hypotheses of 
the theorem. But it will provide an insight as to what properties of the 
functional are to be considered and how to modify the theorem to suit 
such functionals. This will be done in the next section. 
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Theorem 2.5.1. Let n and v be integers > 1. Let Q c R" be an open 
subset and let aeR. (If ,meas Q, =+ oo, assume a = 0). Let 

g : QxR v ^ [a, +oo] 

be a mapping such that 

g(x, .) : R v -» [a, +oo] 
w convex a«<i continuous for almost all xeQ.;, 

g(.,q) : Q -» [a, +oo] 
w measurable for all qeW. Then, if qi — > g weakly L ! (Q) 

(2.5-1) J" g(x, q(x))dx < lim inf J* g(x,qk(x))dx. 

n n 

/« o?/ier words, f/ie mapping 

q -> J" g(x))<ix 
a 

w weafcfy Zower sera/ - continuous on L 1 (O). 

/Voo/ First of all, without loss generality, it can be assumed that g > 0. 
(If meas Q < +oo, replace g by g - a meas (Q)). The continuity in g 
for almost all x and the measurability in x for all g implies that g is a 
Caratheodory function and so if q(x) is measurable in x, so is g(x, q(x)). 
Since now g > 0, the integral 



g(x, q(x))dx 

a 

makes sense. 

Let q n — > g in L ! (Q) strongly. Let [q„ k } be any subsequence such 
that the sequence 

J q(x, q„ k (x))dx 
n 
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is convergent. Now there exists a further subsequence (again denoted 
by q„ k for convenience) such tha q nk {x) — * q(x) for almost all x. Hence 

g(x,q nk (x)) -> g(x,q(x)) 

for almost all x. Thus by Fatou's lemma 

^ g(x, q(x))dx < lim inf ^ g(x, q„ k (x))dx 
n n 

= lim I g(x, q nk (x))dx. 
p. 

Then as the subsequence {q„ k } was chosen arbitrarily subject to the 
condition that the integrals converge the relation ( 12.5-lt follows, when 
<7„ — > q strongly in L (Q). Thus the functional 

(2.5-2) J(q) d = J" g(x,q(x))dx 

Q 

is strongly lower semi - continuous. It is also easy to see that J is convex. 
Now if aeR, then 

{qeh 1 (Q)\J(c) < a] 

is strongly closed and convex and hence weakly closed {Mazur's Theo- 
rem). Thus it follows that J is weakly lower semi - continuous which is 
equivalent to ( I2.5-1I ). □ 

Remark 2.5.1. If g is independent of x, then it suffices to assume that g 
is convex and continuous. 

Remark 2.5.2. The above result conuld be applied as follows: If g(x, q) 
> c + b\q\ p , b > 0, p > 1 then a minimizing sequence will be bounded in 
L P (Q). Since p > 1, a weakly convergent subsequence in L P (Q) can be 
extracted. If Q is bounded, this implies weak convergence in h\Q.) and 
an application of the above result would show that at the limit, / attains 
its minimum. 
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Remark 2.5.3. If g did not take the value +oo, the convexity in q also 
implies continuity. However, if g assumed the value +oo continuity no 
longer follows from convexity. The inclusion of the value +oo in the 
range of g is necessary for applications. It will be needed (as was men- 
tioned in Section l2~T1 i that the energy tends to +oo as det (F) approches 
through positive values. 

Theorem 2.5.2. Let Q. c R" be a bounded open set. Let #:flxR v -> 
[a, +oo] be such that 



for all qeW and almost all xeQ. Let I : 'W ,p {Q) ->lfeci continuous 
linear functional. Let Yq cT be of strictly positive da - measure and let 
U be a weakly closed subset of 



<W{x, .) : R v ^ [a, +oo] 



is convex and continuous for almost all xeQ., 



'Wi^q) : n [a,+oo] 



is measurable for all qeW. Let there exist c, b, p such that 



(2.5-3) 



b > 0, p > 1, and W(x, q) > c + b\q\ p 



(2.5-4) 



V - [veW hp (n)\v = 0on T ) 



Define 



(2.5-5) 





(2.5-6) 



I{u) = MI(v) 



has atleast one solution. 



88 



2. Some Mathematical Aspects of Three-Dimensional . 



Proof. Let v k be a minimizing sequence in U, i.e.v^eU and 
I(y k ) -> inf I(y) < +00. 

Now, 

/(v) > Cmeas(n) + b\v\[ pCl - M\\Mh, P # 
> Cmeas(Q) + ft'|v|f pn - |K||||v||i, p , n 

where b' > 0, using Poincare's inequality. Since p > 1, it follows that 
/(v) — > +00 as ||v||i iPi n +°°- Hence 

sup ||v*||i )Pi n < +oo- 
k 

As p > 1, V ' p (0) is reflexive and a weakly convergent subsquence 
can be extracted. Denoting this subsequence again by v k , let v k — > m 
weakly in W 1,p (£2). Since U is weakly closed, weU. 

Now clearly £(y k ) -> ^(w). Also Vv^ -> Vw weakly in L^(Q) and 
hence (Q. is bounded, in L ! (n). Then by Theorem 12.5. II it follows that 

inf /(v) < /(«) < liminf/(v*) - inf /(v). 

veU i— »oo veU 

□ 

Remark 2.5.4. The application of this result to the linearized elasticity 
system is easy (cf. Exercise 2.5-1. But unfortunately, it is not directly 
applicable to non - linear elasticity. 

The assumptions to be satisfied are firstly, p > 1, which leads to the 
choice of the appropriate space G l p (n), and, secondly, the convexity of 
the function < W. 

If the fuctional / were strictly convex, this would imply uniquer- 
ess of solutions which is not physically acceptable (cf. Section 12. U . 
However, it is not even possible to have F — > 'W(F) convex in three - 
demensional elasticity (cf. Exercise 2.5-4). Note that in the linearized 
elasticity system, it is true that "W is convex, but then one can show (cf. 
Exerise 2.2-1) that a linear model also contradicts the axiom of material 
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frame indifference. Finally, to show U to be weakly closed in W i,p (Q.), 
it is usually shown that it is strongly closed and convex. But typically, U 
will be non - conver with constraints like det (Vt/0 > or det (Vi/r) = 1 
Thus to overcome these two difficulties, the notions of polyconvexity 
and compactness by compensation will be introduced in the next section. 

Exercises 

2.5-1 Show that in the linearized system of elasticity, unilateral condi- 
tions can also be taken into account (apply Theorem 12. 5 . 21 with 

U = {veY\u 3 > on r[, c T - r }, 

where 

V = {veH 1 (Q)|v = on F' }, da meas T > 0). 
2.5-2 For a St Venant-Kirchoff material, 

W(V V ) = -HE) 2 + p tr(E 2 ), A >0,p>0 

where 

E = E(v) = ^(Vv r + Vv + Vv r Vv). 

Show that the corresponding energy is coercive on the space 
V = {veW M (Q); v = on T ] , da -meas r > 0. 

1 T 

2.5-3 If E = -(F lT F - I), show in the above case that F -> W(F) is 107 
convex. 

2.5-4 Show that the convexity of the function F — > "W{F) is physically 
unrealistic (cf. TRUESDELL and NOLL [1955]). 

2.5-5 For a St Venant-Kirchhoffmaterial, show that the solution u ob- 
tained via the implicit function theorem minimizes locally the en- 
ergy in W loo (Q) but not necessarily in W 14 (Q). 
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2.6 J. BALL'S Polyconvexity and Existence Theo- 
rems in Three Dimensional Elasticity 

In the last section, it was seen that the lack of convexity of the stored 
energy function was an obstacle to the application of the existence the- 
orem (cf. Theorem 12. 5 .21 . Now, an extension of the notion of convexity 
following J. BALL will be introduced. 

Recall that if A is a matrix, then adj(A) stands for the transpose of 
the matrix of cof actors of A. The following identity holds: 

(2.6-1) A(adj A) = (adj A)A = det(A)/. 

Thus, if A is invertible 

(2.6-2) adj (A) = det(A)A _1 . 

Also, 

(2.6-3) adj(Afi) = adj(fi) adj(A). 
and 

(2.6-4) (adjA) T = adj(A r ) 

In the study of deformations (cf. Section II. 1L it was seen that 
lengths were modified by a function of F(= V0) via C = F T F. Sur- 
108 face areas were changed in terms of adj(F) and volume elements were 
altered by a factor of det(F). Since it is natural to expect that a stored 
energy function somehow takes these into account, it is reasonable to 
assume that 

(2.6-5) "W(F) = §{F, adj(F), det(F)) 

for all FeM^., where 

§ : M'xM^x]0,ooh R 

is a given function (since for FeM^., adj(F)£M^). While it not true that 
*W as a function of F is convex (cf. Exercise 2.5-4), it is no longer 



2.6. J. BALL'S Polyconvexity and Existence Theorems.. 



91 



impossible to expect Q to be a convex function of its three arguments, 
F, adj(F) and det F (of course, the domain of definition of Q is not a con- 
vex set, but that is easily handled as will be seen below). For instance, 
*W(F) - detF is not a convex function but Q{8) - 8 is convex! 

Let V be a vector space and U c V ary subset. Let J : U — > R be 
a function. It is said to be convex if there exists / : co(U) — > R, where 
co(U) is the convex hull of U, such that 

(2.6-6) /(v) = /(v) 

for every vet/. 

Let U be a non - empty suset of M 3 and let 

(2.6-7) nA = {(F, adj(F), detF)\FeU\. 

Thus <U c M 3 x M 3 x R. A function TV : 'W R is said to 
be polyconvex if there exists a convex function Q : H — > R such that 
( !2.6-5t holds for every Fell 

If U = M 3 then co(U) - M 3 and U = M 3 x M 3 x]o, +oo[ while co 
(^) = M 3 x M 3 x]o, +oo[. Thus a stored energy function *W : M 3 R 
is polyconvex if there exists a convex function 

: M 3 x M 3 x]o, +oo[— > R 

such that d233l holds for all FeM 3 . 

The condition of polyconvexity on the stored energy function leads 
to a class of hyperelastic materials known as OGDEN'S materials 
(which we now define for the compressible case: for the incompress- 
ible case, see Exercise 2.6-8). 

As a simplest possible example, let a > o, b > o and T :]o, +oo[— > R 
be a convex function. Define 

(2.6-8) <W{F) = a\\F\\ 2 + b\\ adj F\\ 2 + T(det F) 

where 



(2.6-9) 



||F|| 2 = tr(F T F) = F:F. 



92 



2. Some Mathematical Aspects of Three-Dimensional . 



Clearly as ||.|| 2 is a convex function, and as Y is also convex, it fol- 
lows that *W denned by ( !2.6-8t is polyconvex with 

(2.6-10) Q{F, H, 6) = a||F|| 2 + b\\H\\ 2 + T(S). 

Remark 2.6.1. It will later be assumed that T(6) — > +oo as 6 — > o + . 

Let FeM+. Let U = (F T F) 1 ^ 2 with eigenvalues vi, v%, V3. These are 
called the principal stretches of T 7 . Then it is easy to see that 

(2.6-11) ||F|| 2 - v\ + v\ + v 2 

(2.6-12) || adj F|| 2 = v\v\ + v 2 v 2 + v\v\. 



Hence (!2.6-8t will now read as 

(2.6-13) <W{F) = a(y\ + v\ + vf) + J>(v?vij + v^v 2 , + v?v 2 ) + r(viv 2 v 3 ). 

This can be generalized to get an OGDEN Material as follows. Let 
T :]o, +00 [— > R be a convex function. Let a, > 0, 1 < i < M, bj > o, 1 < 
j < N. Further let, 



(2.6-14) 

Now, if C - F T F, 
(2.6-15) 



1 < a\ < ■ ■ ■ < ccm 

1 <y6l <---<P N . 



(tr(C a/2 ) = v a x + v a 2 + 

(tr((adj Cf> 2 ) = (viv 2 f + (v 2 v 3 y ? + {v z v { f. 



Now define for FeM;;, 

M N 

(2.6-16) W(F) = J] m tr(C a ' /2 ) + ^ tr(adj C/*' 72 + T(det(F)). 

i=l 7=1 

Theorem 2.6.1. (/J An Ogden's Material is polyconvex. 
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(ii) It satisfies the following coerciveness inequality. 

(2.6-17) <W(F) > C + Ci\\F\\ aM + C 2 || adj(F)f" + T(det(F)). 

Proof. Each summand in d2.6-16t is a symmetric function of the eigen- 
values v,t(resp. VkVk+i) and is convex and non-decreasing with respect 
to each variable on ]o, +oo[ 3 . Such a function is convex with respect to 
F (resp. adj F) and thus 'W is polyconvex. (cf. Exercise 2.6-1). □ 

To prove the coerciveness notice that 

tr(C a/2 ) = v a 1 +v a 2 +v a 3 

> C(ar)((v 2 + v\ + v 2 ) 1/2 f 
= C(a)||F|| a 

for any a > 1, by the equivalence of norms in R 3 . The inequality 
(I2.6-I7t follows directly from this observation. 

The stored energy function of a St Venant-Kirchhoff material: 

(3A + 2u\ U + 2u\ , /I 

(2.6-18) W) = -(—4-^) W + I" 8"^J tr(C } + 4 tr ( ad j(0) 

(6n + 9A\ 



is not polyconvex (Exercise 2.6-2). This stems from the fact that the 
coefficient of tr{C) is < 0. 

Remark 2.6.2. In the incompressible case, det(F) = 1. Thus *W(F) 
has no dependence on det(F). Such Ogden's materials comprise the 
so-called MOONEY-RIVLIN materials. 

To fix ideas, consider an Ogden's material described by 

M N 

(2.6-19) *W(F) = J] a t tr(C a ' /2 ) + J] fe/(/r(adj Cf" 2 ) + T(det F) 
with a\ > o,l < i < M,bj > o, \ < j < N and T(5) > C5 r + J, c > o. 
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Then by Theorem l2.6.11 the following coerciveness condition holds: 

(2.6-20) W(F) >a + b(\\F\\ p + || adj F\\ q + (det F) r ), p 

= max at,q- max XyS/. 

112 Since it is natural to desire that 

r W(F)dx < +oo, 



r 



the set of admissible deformations will typically be of the form 

(2.6-21) U - jiAeW 1 ' p (0)| adj (VtA) eh q (Q), det (VtA) £Z/(Q) det 

(Vi/r) > o a.e., if/ - (p Q on T } . 

Such sets are 'highly' non-convex (cf. Exercise 2.6-3). It will not 
be shown that U is weakly closed (in fact that is not true in general) but 
it will be shown that weakly convergent subsequences can be extracted 
which suit the purpose of minimizing the energy. 

In order to do this, the mappings (p — > adj(V0) and cf> — > det(V^) 
have to be looked at more closely. 

Counting indices modulo 3, the matrix adj(V0) can be defined by 

(2.6-22) (adj (V0)) i; - = {d M ^j^id i+2 <p j+2 - di +2 ^j^id M <p j+2 ). 

If (t>eW Up (Q),p > 2, then it is easy to see that adj(V0)eL p/2 (Q). 
The mapping defined in this fashion between W 1,p (0) and L p ^ 2 (0) is 
non-linear and continuous. We denote weak convergence by — s 

Theorem 2.6.2. Ifcf>eW UP (Q.),p > 2, then adj(V0)eL p/2 (Q) Further 

6 n d>inW hP (n),p > 2 ) 
(2.6-23) } implies H = adj(V0). 

adj(V0") H in L g (Q), q > 1 j 

113 Proof. (1) First an alternative definition of adj(V0) will be established 
in the sense of distributions. Let (f>eC°°(Q.). Then a simple computation 
yields 

(2.6-24) (adj(V0)) y = d i+2 {<p j+2 d M <p M ) - d M (<f>j +2 d i+2 <p j+l ) 
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(with no summation on i and j). If 9eD(Q.), then 



(2.6-25) f (adi(V cf>))ij9dx = - [ 




+ (f>j +2 dj + 2(f)j+idi + i9dx 
Jn. 



For fixed 9eD{Q), a simple application of Holder's inequality shows 
that each term in (!2.6-25t is a continuous function if 0eH (Q) (where if 
4>eW i ' P (Q.), p > 2). Since C°°(Cl) is dense in each of these spaces, it 
follows that (12.6-251 is true for 4>eW UP (Q.),p > 2. Thus (12.6-241 holds 
for 0eW 1,p (O) in the sense of distributions. 

(ii) Let <p n cf> in W UP (Q). Let 9 = (9^,9^(0) . Let p* = +oo 
if p > 3 and be given by 



for p < 3. Then for 1 < q < p*, 

(2.6-27) W 1 ' P (Q) ^ L q (Q), 

i.e., the above inclusion is compact. Now if^eD(Q) is fixed and ifr" — 1 i/f 
in W 1 " P (Q) and <p n -± <f> in W^O), then 



i/r" ^ i/r in L^Q) 
-> d k ^>mL q {Q). 



If further — h 
P 



— < 1, it will then follow that 
1 




From this observation and from (12.6-25 l l. it follows that 



(2.6-28) 



f (adj (V<f)) : &/* -> f (adj (V0")) : to 
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provided that q < p* and — I — < 1. It is easy to see that this is 

q p 

equivalent to p > 3/2, which is satisfied anyway. 

(iii) Let (f> n ->■ <f> in W hP (Q) and adj(V0") H in L 9 (H). By (ii) 
above for any = (0,y), e£)(il), (12.6-281) holds and also 

| (adj (V<f)) : 6dx -» f // : 0dx. 

Thus 

(2.6-29) f (adj (V0 n )) -H):6-o 

Jo. 

for all 6» - (0 y ), 0, 7 e£>(Q) and since (adj(V0")) - //)eL 1 (Q) and D(Q) is 
dense in L ! (Q), it follows that 

H = adj (V0) 

thus proving the theorem. □ 

Remark 2.6.3. If p* is as in H2.6-26t . if <j>jeLP*{Ci) and di<p k eLP(Q.), and 

if — I < 1 (which is equivalent to p > 3/2). the product (pidjcpk 

P P* 
belongs to L l (Q.). 

Now using (!2.6-25b . the adjugate of Vcp for (peW l,p (Q.) can be de- 
fined in the sense of distributions. This definition of the adjugate of V<p, 
extended to p > 3/2 is denoted by 

Adj (V0) . 

The step (ii) of the proof of the above theorem goes through for 
p > 3/2 as remarked in the proof itself. Hence for p > 3/2, if <f>" — ^ <p 
in W liP (Q), then Adj{V<f> n ) (V<p) in the sense of distributions. 

Remark 2.6.4. The above theorem implies that the set 



U = {((f>,H)eW l ' p (0) x L q (0)\H = adj (V0)) 
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is weakly closed in W hP (Q.) x U'(Q) for p > 2, q > 1. But the set 

{4>eW UP (pi)\ adj (V0) eh q (Q.)] 
is not necessarily weakly closed in W l ' P (Q.). (cf. Exercise 2.6-4). 

Now consider the mapping cf> — > det(V0). In the first place 
(2.6-30) det(V0)=^ sgn ^d^^d^^d^aO) 

<T£P 3 

where is the set of all permutations of (1,2,3)- Now det(V^) will 
be in L l (Q) if (f>eW 1 ' P (Q.), p > 3. But one use the integrability of the 
adjugate and improve on this by noting that 

(2.6-31) det(V0) = 5 I -0i(adj(V0)) il 

(summation with respect to i). Now, if drfi e L P (Q.) and (adj(V^)),i e 
W(Ci) and if l/p + l/q < 1, then det(V^)6L 1 (^)- Thus (I2.6-31H will 
be used to define det(V0), for if p < 3, the formula (l2.6-30t makes no 
sense. 

Theorem 2.6.3. Let p>2, <f>eW UP (Q) such that adj(V^)eL p '(Q), where 

— I = 1. Then det(V0) given by (12.6-311) is in L 1 (Q). Further 

P P' 



(2.6-32) 

<p" (f>inW l ' P {Q.),p > 2 

adj (V0 n ) H in h q {Q), - + - < 1 

P q 

det(V0") — * 8 in L r (Q), r > 1 



implies 



\H = adj (V^) 
[5 = det(V0). 



Proof, (i) The main difficulty in the proof is to give an alternative def- 
inition of the determinant in the sense of distributions. Let 0eC°°(n). 
Then 



det(V0) - <3i0i(adj(V0)) n 
= d,(0i(adj(V0)) a ) 
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using the Piola identity (cf. Exercise 1.1-1). So, for 9eD{Q) 

(2.6-33) I det(V0)ft/jc = - ( 0i(adj(V0)) n d/0dx. 

Jn Jn 

It will be shown that (12.6-331) is valid for (p e W 1 '^) with adj(V0) 

» 1 1 
€ 1/ , where — I 1. 

Now by the Piola identity, for smooth 0, <9;(adj(V0))n = o or for 
0e£>(Q) 

(26 ' 34) JwW»aW* = o. 

Q 

By the density of smooth functions in W 1,P (Q), it follows that 
(12.6-341) is true for all <f>eW UP (n),p > 2. 

For expository convenience, set <f> - <j>\, W,- = (adj((V0)),i. Then 

0ew 1 ' P (D) and w/eL^n), - + — = 1. 

P P 

Let pe£>(R 3 ),p > o and f p = I. Define 



(2.6-35) p k (x) - £ 3 p(£*) 

so that pk has the same properties as p together with the property that 
supp (pk) skrinks to zero as k — > oo. Let w* be extended by o outside f2 
and define 



(2.6-36) {p k * Wi){x) = | p(k)(x - y)wi(y)dy. 

Then the function p* * w\ is smooth and converges to w, in L p (Q.) 

Let 6eD(Q.) be fixed. Then there exists & = k o {0) such that the 
support of the map y — > p^(x - y) is contained in O for all k > k and for 
all ;cesupp(#). then for k > k , 

r d 

(2.6-37) div(p* * w)(X) - I -7j-pk(x ~ y)wi(y)dy = o 

n 
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using A2.6-34I I. for xfsupp(6>). Now if /eC°°(Q) and / -> (p in 
W lP (Q.), it follows that 

- J 4> k {pk * Wi)dj6dx - J di(f> k (pk * Wi)6dx + J 4> h di{pk * Wj)9dx. 

fi supp(ff) 

By d2.6-37t the second integral on the right-hand side vanishes. 
Now passing to the limit in each of the integrals as k — > oo, 

- f <pWidi9dx = ^ di4>Wi9dx. 
n n 

from which <!2.6-33b follows. 

(ii) If (p n -> <p in W 1 ' P (Q) and adj(V0") — adj(V<£)in L*(fi), with 

p > 2 and — + - < 1, then using the same type of compactness argument 
P 1 

as in the proof of theorem r2.6.2l it can be shown that for 9eD(Q.), 

J det (V0„) 6dx -> J det (V0) Qdx. 
q. n 

(ih) Using the previous step, the conclusions of the theorem can be 
drawn exactly as in Theorem l2.6.2l □ 

Remark 2.6.5. As in the case of adj(V0), if p > 3/2 the determinant can 
also be defined in £>'(Q) using the fact that (f>i(adj(V4>))n€L 1 (D.). The 
distribution obtained is denoted by Det (V0). 

Remark 2.6.6. The above theorem shows that the set 

{(</>, H, 6)eW hP (Q) x L q (Q) x L r (Q)\H = adj (V0) , 5 = det (V0)) 

is weakly closed in / '•''(£}) x U(Q). But the set 

{0eW 1 " p (Q)| adj (V0) eL 9 (Q), det (V0) eZ/(Q)j 

is not necessarily weakly closed in W l ' P (Q.). 
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The following existence theorem can now be proved. 

Theorem 2.6.4 (J. BALL). Let <W : M 3 -> R be a stored energy func- 
tion, such that 

(i) (Polyconvexity) there exists Q : M 3 x M 3 x]o, +oo[— > R which is 
convex and such that for all FeM 3 + , 

(2.6-38) <W{F) = §{F, adj(F), det(F)); 

( ii) ( Continuity at +oo ifF n — > F in M 3 , H n —> H in M 3 and 8 n — > o + , 
then 

(2.6-39) limg(F n ,H n ,6 n ) = +oo; 



(Hi) (Conerciveness) There exist aeR,b > o,p > 2,qeR with — h - < 

p q 

1, and r > 1, and r > 1, such that and for all {F,H,S)eM? x 
M 3 x]o, +oo[, 

(2.6-40) Q(F, H,8)>a + b(\\F\\ p + W\\ q + S r ). 

Let QcR 3 kfl bounded open subset with boundary T = T U T\ 
where the da-measure ofT a is > 0. 

Let /eI/'(Q),,£>e]L' T (ri) such that the maps 

(AeW 1 " p (Q) J f ■ ipdx and if,eW l ' P (Q) J g.tyda 



are continuous. 

ThP, 



Let I : W l *(Cl) ^Rbe given by 
(2.6-41) I(\fj) = J W(Vi[,)dx-\J f.ifjdx + J g.ipda 



n vn 
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and U c W l ' p (0)be defined by 

(2.6-42) U = j(A£W 1 ' p (Q)| adj (V<f>) eL q (Q), det (V^) e7/(Q), det 

(VtA) > o a.e., tA = <p on Y (> ) , 

with o eW 1,p (Q). Assume that U - (f> and that 

inf 7(i/0 < +oo 

l/reU 

TTjera ?/ze problem: Find fieU such that 
(2.6-43) 7(0) = inf 70/r) 

\fie\S 

has at least one solution. 

Proof. Step (i). Transformation of the problem. Define § : M 3 X M 3 X 
R — > RJ/j+oo) by 



(2.6-44) Q{F,H,S) = 



\§{F,H,8)\fS> o 
+oo if S < o. 



Then Q is easily seen to be convex and continuous into [a, +oo]. 
Thus the functional 

(2.6-45) I(ifr) = J Q (V^, adj (VtA) , det (Vi/0) 
n 

- J > f.if/dx + J" g.i//da 
,n ri 

is well-defined. 

Let if/eW^iQ) with adj(V^)eL«(Q) and det(V^)eL r (fl). If 7(0) < 
oo, then it follows that det(VtA) > o a.e. 

If iJ/eU then 7(0) = 7(0). Thus the original problem is equivalent to 
minimizing I(if/) over U. 
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Step (ii). It can be shown that (cf. Exercise 2.6-5) for all i//eW lP (Q) 
with ip = <p on r o , 



(2.6-46) 
where d > 0. Hence for i//elJ, 



J \if/\ p dx<d 

Q 







P 


f \Viff\ p dx + 


r \(po\da 




Si. 







I(i/f) > a meas Q + b' J + b J || adj (V^) \\ q dx 

+ 6 J (detWy-CH.Alli^n 

Q 

with b' > o, or 



(2.6-47) /OA) > Co + C 1 || t A||^ n + C 2 |adj(V ( A)|^ n + C 3 |det(V t A)i; rn 

121 with Ci , C 2 , C 3 > 0. (cf. Remark l2~6~7l) 

Step (iii). Let <p n € U be a minimizing sequence for /. From the coer- 
civeness (!2.6-47t . it follows that <p n , adj(V0") and det(V0") are bounded 
in W l ' p (Q.), L 9 (Q) and Z/(Q) respectively. Since these spaces are reflex- 
ive, a subsequence <p n can be found such that 

(f>" -» (f> W ! ' p (0) 
adj(V0") H in L 9 (Q) 
det(V0") 8 in U(Q). 

But by the previous theorem H = adj(V0) and 6 = det(V0). Thus 
by the convexity and contnuity of Q, I is weakly lower semi-continuous 
(cf. Theorem l2.5.U and so 

7(0) < liminf/(0") < +oo. 

n— > oo 

Hence det(V0) > o. It can be show that 4>\r„ = <Po (cf- Exercise 
2.6-6) and so (f> e U and it follows that I and hence / attains a minimum 
at (p. □ 
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Remark 2.6.7. The coerciveness condition ( 12.6-471 can be obtained 
only on U and not on the whole space. It uses the fact that <p = <p 
on r o for all the functions under consideration. 

Several comments on the above result are in order here. First of all, 
unlike the approach based on the implicit function theorem, the result 122 
is applicable to "all" forces (not just "small" ones) and to all boundary 
conditions. Of course, in case of the pure traction problem, the forces 
must satisfy certain compatibility conditions. It is also applicable to the 
mixed displacement-pressure problem (cf. Exercise 2.6-7). 

A shortcoming of this approach is the lack of regularity of the solu- 
tion. Here it is not know if the minimizing function satisfies the equilib- 
rium equations even in a weak sense. Further even thought it is true that 
the solution satisfies det(V^) > a. e. , additional conditios are needed 
to insure that <p is one-one (see BALL [1981c]). 

It is possible to extend this approach to cover the incompressible 
case where det(V^) = 1. (cf. Exercise 2.6-8). 

Consider a St Venant-Kirchhoff material. If the forces are 'small 
enough' it was shown that there exists a 'small' solution to the pure dis- 
placement problem (Theorem l2.3.H . However, in the pure displacement 
or mixed displacement traction problem, owing to the non-polyconve- 
xity, it cannot be shown that the energy is minimized. An open prob- 
lem is to prove existence of 'small' solutions for small forces the mixed 
problem such materials. 

To conclude this section, it will now be examined how to choose 
a stored energy function given a compressible material. Consider a 
compressible material (steel, for instance!). Around a natural state it 
is known that the stress tensor can be written as 

(2.6-48) S* (E) = A{tr E)I + 2pE + o(E), 

where E is the Green-St Venant strain tensor. The constants A and p are 123 
strictly positive and can be determined form experiments, albeit approx- 
imately. It has been shown that (cf. Section [L4l 



(2.6-49) 



<W(F) = -(tr(£)) 2 + ptr(E 2 ) + o(\E\ 2 ). 
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The foillowing theorem due to CIARLET and GEYMONAT [1982] 
says that it is possible to express such a material as a simple Ogden's 
material. 

Theorem 2.6.5. Given A > 0, p > 0, it is possible to find a > 0, b > 
and a function T :]0, +oo[— > R which is convex satisfying 

(2.6-50) Y(8) >C3 2 + d,C>o 

such that the corresponding stored energy function 

(2.6-5 1) 'WiF) = a\\F\\ 2 + b\\ adj F\\ 2 + T(det F) 

agrees to /l/2(tr(£)) 2 + p tr(£ 2 ) upto o(\E\ 2 ). 

Proof (Sketch). Setting C = F T F = I + 2E, then 

||F|| 2 = trC - tr(/ + 2£) 
||adjF|| 2 = tr(adj(/ + 2£)) 
det(F) - V( det (/ + 2E)). 

Expanding these about I, it it follows that 

W(F) = 3a + 3b + T(l) + (2a + 4b + T'(l)) tr(£) - (2b + T'(l)) tr(£ 2 ) 
+ (2b + i(T(l) + T"(l)))(tr(£)) 2 + o(\E\ 2 ). 

Comparing with d2.6-49t . it follows that 

(2.6-52) 2c? + 46 + r'l=0, 

(2.6-53) -(2b + T'(l)) - p. 

(2.6-54) 26 + i(T'(l) + r(l)) = |. 

These equation must be solved such that a > o,b > o and T"(l) < o 
(T is convex). (By ( 12.6-521 )) it follows that F(l) < o). It is easy to 
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see that any point (T'(l)), r"(l)) on the open line-segment shown in 
Fig. l2.6.l1 give a > o,b > o satisfying the above equations. 



r"(i) 




Figure 2.6.1: 

Now to choose a convex function T :]0, +oo[— > R. satisfying d2.6-50t . 
One can find a > 0, /? > such that 

(2.6-55) T(6) = aS 2 -y31og<5. 

This function also is such that Y(6) — > +oo as 5 — > o + . □ 

It follows now that the associated minimization problem has at least 
one solution by J. BALL'S theorem. Here 

(2.6-56) U - {ifr e H'(Q)| adj(VtA) e L 2 (Q), det(VtA) e L 2 (Q), det(V^) 

> Oa.e. and ij/ = <po on To} 

Remark 2.6.8. In ( I2.6-50I) . the term C<5 could have been replaced by 125 
C5 r , for r > 1 . The definition fo U would be modified accordingly. 

Remark 2.6.9. It is also possible to choose X{F) in the form 

(2.6-57) JfT(F) = a x \\F\\ 2 + a 2 \\F\\ 4 + b\\ adj F\\ 2 + T(det F) 

where a\ > 0, ci2 > 0, b > and T convex. (The St Venant-Kirchhoff 
stored energy function resembles this, only a\ < 0.). In this case, it can 
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be seen that the admissible range of values (r'(l), r"(l)) lies in the open 
triangle of Fig. l2~6.2l (Exercise 2.6-9). 



r"(i) 




Figure 2.6.2: 



Exercises 



2.6-1 Let Jtf : M 3 -> R be a function such that 

<H(F) = (f>(vi,v 2 ,v 3 ),F € M 3 + 

where Vj, i = 1, 2, 3 are the principal stretches of F. If (p is a sy- 
metric function which is convex on (]o, +oo[) 3 and non-decreasing 
in each variable, show that 'W is convex. 

2.6-2 Show that the stored energy function *W for a St Venant-Kirchhoff 
matrial (cf. ( 12.6- 18t ) is not polyconvex. 

2.6-3 Show that the set 

e W 1,p (Q)| adj(V^) € V{Q),q> 1} 

is not convex (p > 2). 

2.6-4 If p > 2 and q < p/2 show that (f> n -» <p in W'^CQ) implies 
adj(V0") adj(V0) inL«(Q). 
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2.6-5 Show that there exists a constant d > such that for all if/ e 
W lp (Q.),i// = tAo onr , 



a 



J* \ij/\ p dx<d J \Vi/f\ p dx + J 



\<f> \da 



2.6-6 If f (f> in W 1>P (Q) and 0" = O on r o , show that cf> = <p on r o . 

2.6-7 Apply J. BALL's theorem to the mixed diaplacement-pressure 
problem. 

2.6-8 . Let U be defined by 

U - {{(f>, H) € H^Q) x h 2 {Q)\H = adj(V0), <p = <p Q on T Q , 

det(V0) = la.e} 

(incompressible case). Assume U = (i) Show that 7/ is weakly 
closed in the product space H^Q) x L 2 (Q). (ii) Consider 



*W{F) = a\\F\\ 2 + b\\ adj F\\ 2 , a>o,b>o, 
a n ri 



Show that the problem: Find <p e U such that 

U = {4> e H 1 (Q);adj V0 e L 2 (Q),0 = 4> ) on r o ,detV0 = la.e.) 

7(0) = inf/W 

has at least one solution, (iii) If is smooth, show that the La- 127 
grange multiplier arising out of equality constraint det(V0) - 1, 
is the pressure, (cf. Exercise 2.1-2). 



2.6-9. Check the assertion made in Remark l2.6.9l 
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128 Bibliography, Comments and some Open Problems 

No attempt has been made to give an exhaustive list of pertinent 
references. 

The first chapter of these lecture notes gave a description of elastic- 
ity in three demensions. For further refrences, one may also consult 
GERMAIN [1972], GREEN and ZERNA [1968], GREEN and AD- 
KINS [1970], GURTIN [1981a, 1981 b], MARSDEN and HUGHES 
[1978,1983], STOKER [1968], TRUESDELL and NOLL [1956], 
VALID [1977], WANG and TRUESDELL [1973], ERINGEN [1962] 
and WASHING [1975]. 

The second chapter discussed some methods for proving the exis- 
tence of solution to the boundary value problem of non-linear elasticity 
and to the associated variational problem, in the case of hyperelastic 
materials. 

For references about the linearized system of elasticity, see DU- 
VAUT and LIONS [1972], FICHERA [1972] and GURTIN [1972]. 

The key result in proving existence via the implicit function theorem 
is the ly^C^-regularity of the linearized system of elasticity. The case 
p - 2 was proved by NECAS [1967] and the regularity for other p was 
proved by GEYMONAT [1965]. From this regularity result, (proved 
however only for the pure displacement problem) the existence theorem 
was independently proved by CIARLET and DESTUYNDER [1979b], 
MARSDEN and HUGHES [1978], VALENT [1979]. The basic idea, 
however, goes back to SPOPPELLI [1954] and VAN BUREN [1968]. 
The extension of this result to more general constutive equations was 
particularly studied by VALENT [1979]. See also VALENT [1978a, 
1978b]. 

129 The necessity of the W 2 p (Q)-regularity of the linearized problem re- 
stricts the application of this method to pure displacement problems. It 
is also possible to treat the pure traction problem, which is more compli- 
cated owing to the compability conditions which the given forces must 
satisfy. For details see CHILLINGWORTH, MARSDEN and WAN 
[1982]. 

The increment method described in Section |2~4"1 is none other than 
Euler's method for approximating an appropriate differential equation 
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in a Solbolev space. In other words, this method appears as an infi- 
nite demensional version of the so called continuation by differentiation 
approach as described, for instance, in RHEINBOLDT [1974]. 

To the best of the author's Knowledge, the convergence of increment 
methods for non-linear elasticity problems has been analysed so far only 
in some special cases, such as the one demensional model of a thin shal- 
low spherical shell, by ANSELONE and MOORE [1966] or some finite 
dimensional structural problems by RHEINBOLDT [1981]. The results 
presented in these lecture can be found in BERNADOU, CIARLET and 
HU [1982]. 

For a decription of incremental methods in non-linear elasticity see 
MASON [1980], ODEN [1972] and WASHIZU [1975]. 

The variation approach is based on the famous article of BALL 
[1977]. In addition to the notion of polycomvexity another essential 
contribution of J. BALL is that one can pass to the weak limit in cer- 
tain non-convex sets as was seen in Section 2.6. This idea of compact- 
ness by compensation was also developed by MURAT [1978, 1979] and 
TARTAR [1979]. See also AUBERT and TAHRAOUI [1982]. 

Other important refrences are BALL [1981a,1981b, 1981c], BALL, 
CURRIE and OLIVER [1981], BALL, KNOPS and MARSDEN [1978]. 
See also EKELAND and TEMAN [1974] for the general problem of 
minimizing functionals. 

The notion of polyconvexity led to the definition of an Ogden 's ma- 
terial (cf. OGDEN [1972]). A St Venant-Kirchhoff material is not Og- 
den's material and the existence of a solution to the corresponding mixed 
displacement-traction problem is open. In this connection see also AT- 
TEIA and DEDIEU [1981] and DACOROGNA [1982a, 1982b]. For yet 
another approach, see ODEN [1979]. 

One of the drawbacks of J. BALL'S approach is the lack of regularity 
of the solution and so one does not know if the solution thus obtained 
satisfies the equilibrium equation even in a weak sense. In this contex, 
see the results of LE TALLEC [1981] and LE TALLEC and ODEN 
[1981] for incompressible materials. 

To conclude, we present a list of some of the open problems in non- 
linear elasticity. Some of them have been mentioned in the text before. 
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1. Let C = V0 V0. If C - / is 'small', in a sense to be made precise, 
can it be said that <p is close to a rigid deformation? If some boundary 
conditions are imposed, can it be shown that tp is one-one? 

In this context cf. KOHN [1982], ALEXANDER and ANTHAN 
[1982], ANTMAN [1979]. 

2. The standard implicit function theorem approach fails for mixed 
problems. Could "hard" implicit function theorem like that of NASH 
and MOSER be used? In case of special domains like a thin plate, the 
singularities are know explicity. Could this be used, and the implicit 
function theorem used only on the "regular" part of the solution? 

3. Study of incremental methods taking into account the finite element 
methods. 

4. An incremental method can be formally written down for the mixed 
problem. If it can be shown to be convergent, this would provide an 
existence theorem for the mixed problem. 

5. The minimization procedure of J BALL does not imply that the so- 
lution is small if the forces are small. How can one "distinguish" 
the expected small solution in this case? (In the case of the pure 
displacement problem, the solution via the implicit function theorem 
does not seem to be a local minimum of the energy in the "right" 
space). 

6. A study plasticity has been taken up by TEMAN and STRANG- 
[1980a, 1980b]. They use the linear theory in the part corresponding 
to elasticity. Can one obtain better results by incorporating the non- 
linear theory, using J. BALL'S approach? 

7. A study of 'non local' constotutive equations. Here the comstitutive 
equation is of the form 




where pk is a mollifier. 
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8. One of the hardest open problems of the study of the evolution prob- 132 
lem which is a non-linear hyperbolic problem. The only available 
results are in the one-demensional case due to DIPERNA [1983]. 

See also HUGHES, KATO and MARSDEN [1976]. 

9. Plate theory. A plate can be thought of as a domain £l e - wx]- e 
, + e [, where to c R 2 is a bounded open set and e> o is a small 
parameter, (cf. Fig.l) 




Figure 1 



By methods of asymptotic expansions, the solution (u e ,cr e ) can be 
formally expanded as 

(u e ,o- e ) = (u°,o- Q ) + (u\o- 1 ) + --- 

where (w°,cr ) satisfies a well-known two-dimensional plate model. 
In the linearized theory CIARLET and DESTUYNDER [1980a], 
CIARLET and KESAVAN [1980], DESTUYNDER [1980] have stu- 
died the problems extensively. One can compare the three dimen- 
sional and two dimensional problems and show that (for example) 

H d - 4J\^ ^ o as e ^ _ 

l|w|Ji,n e 



The problem is to numerically verify this. Computing by the finite 
element method, one gets u^ d h and u\ d h approximating u^ d and u^ d 
respectively. Since e is small, unless h is of the same order, the linear 133 
systems become very ill-conditioned. But if h is of the same order 
of 6, the solution u^ d h is not very accurate. Thus to find a better 
method of approximating these solutions and verify the convergence 
described above. 
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10. In the nonlinear case CIARLET [1980] (see also CIARLET and- 
DESTUYNDER [1979], CIARLET and RABIER [1980]) has shown 
that with certain boundary conditions the three dimensional plate 
model for a St Venant -Kirchhoff material is approximated (formally) 
by the well-known two-dimensional von Karman model. While the 
latter has a satisfactory existence theorey, the former has none. If at 
least for e small enough it can be shown that the three dimensional 
problem has a solution converging to a given solution of the dimen- 
sional problem, an existance theorem for such special domains can 
be obtained. 
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List of Notations 



General Conventions: (1) Unless otherwise indicated, Latin indices 
take their values in the set {1, 2, 3}, and the repeated index convention 
for summation is systematically used in conjunction with this rule. 

(2) If a quantity is denoted X in the deformed configuration, the 
corresponding quantity in the reference configuration is denoted Xr. 

Vectors and Matrices 



: orthonormal basis in R 3 
: vector v with components v; 
: matrix A with elements A ;J - (i : row index, 

j : column index) 
: Euclidean inner product 
: Euclidean vector norm 

+ 1 if (i, j, k) is an even permutation of 

(1,2,3) 

- 1 if (i, j, k) is an odd permutation of 

(1,2,3) 

otherwise 
uAv = SijkUjbkei : cross product in R 3 
: matrix inner product 
matrix norm associated with the matrix 
inner product 



(ed 
v = (v,-) 
A - (Aij) 

U - V = U{Vi 

\u\ = V M ' u 



A:B = AijBij = ti(AB T ) 



\\a\\ - ^IaTa 



(A [ y (A . inverse matrix; 
A T : transposed matrix). 



adj A 

l A = (h(Ay,l 2 (A),h(A)) 

h(A) 
h(A) 



adjugate of a matrix (transpose of the 
cofactor matrix) 

set of the principal invariants of a matrix of 
order 3 
an = tr(A) 

j(auajj - aijaij)(- det AtrA -1 if A is 
invertible) 
Z 3 (A) = detA 
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M 3 : set of all matrices of order 3 
M 3 + = {F e idel /•' > 0} 
Q 3 = {F e M?\F T F = FF T = /} 
G 3 = O 3 nM 3 = {F e G 3 |detF - 1} 
S 3 = {F e M 3 |F - F T ) 
S 3 = {F e S 3 |F is positive definite} 
F = RU = VR : polar factorization of an invertible 
matrix (R e O 3 ; U,V e S 3 ) 
C l/2 : square root of a matrix C e S 3 

142 

Functions and Function Spaces 



Id : identity mapping 
v'(a) : Frechet derivative of the mapping v at the 
point a 

Q\a\ v 

(multi-index notation for partial 
derivatives) 
dW idW \ 

~dF (F) = [ dF~ {F) ) 6 M (f ° r a mapping 

W : M 3 ^ R) 
X <—* Y : the canonical injection from X into Y is 

continuous 

c 

X^Y : the canonical injection from X into Y is 
compact 
— 1 : weak convergence 
C°(X, Y) : set of all continuous mappings from X 
into Y 

C m (X; Y) : space of all m times continuously 

differentiable mappings from X into 
F(l < m < oo) 
C m (X) = C m (X; R), < to < oo. 
W m ' p (a) = {v € V(Q); d a v e LP(fl) for all \a\ < to} 
= LP(n)) 
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H m {Q) = W m ' 2 (Q.) 
h p (Q), W m ' p (Q.), H m (Q) : corresponding spaces of vector- valued, or 

matrix-valued, functions 
Mo,p,n : norm of the space L p (fl), 1 < p < oo. 
Hv||m^!o = {/ I \d a v\Pdx} l/ P: norm of the space 

W m ' p (Q), 1 < p < oo 
l|v|| m ,oo,n = max |<9 a v| ,oo,n: norm of the space 

\a\<m 

W m -°°(Q) 

|vk^ = {/ Z |d<*v|^;t} 1/p , 1 < P < oo 

a H= m 
|v| m ,oo,n - max |<9 a v| ,oo,n 

|a|=m 

= (ve c rf°°(Q.y, supp v is a compact subset of 
W™' P (Q) : closure of 9(Q) in 

143 

Miscellaneous 

[a, +oo] - [a, +oo[u{+oo}, d€fi 

f{x) = 0(x) : hm ^ = 
*->o |W| 

c o U : convex hull of U (smallest convex set 
containing U) 

144 Notations in the Deformed Configuration 

S$ = (P(£$r) : deformed configuration 
X = <P{Xr) : generic point of 38 
dSB : boundary of SB 
dS3 = d& Q U d3§\ : JA-measurable partition of d£8 
n : unit outer normal along dSS 
dX : volume element in SS 
dA : surface element on dSB 
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GRAD0 = |^| ( for n 1 ping : 



dT, 



DIV T = -rrrei € ^ 3 : divergence of a tensor field T : 



dXj 

d(X) e R 
6(X> e R 3 
fi(X) e R 3 



f(X,n),X€ 



N = i 



density per unit mass at X e SS 

body force density per unit mass at X € SB 

applied surface force density per unit area 

ofd^atXed^ 

Cauchy stress vector in SB 

Cauchy stress tensor at X e SB 

response function for T = T(F) = T(B), 

with F e M 3 , B - FF T . 



Notations in the Reference Configuration 

SBr, or Q : reference configuration 
Xr, or x : generic point of SBr 
OSBr, or T : boundary of SBr 
OSBr = dSBoR U dS§\R, or T = To U Y\ : <fA#-measurable partition of 

dSB R 

tir, or v : unit outer normal along 8SBr 
cIXr, or dx : volume element in SBr 
cIAr, or da : surface element on 3SBr 
4>, if/ : SBr — > R 3 : deformation of SBr (smooth 

maps with det • > 0) 
u, v : SBr — > R 3 : displacement 

(<p = Id +u, if/ = Id +v) 
d 



DWrTr = -ffi^-ei <= IK 3 : divergence of a tensor field 



7 



dX Ri 
Tr : Y, R 



= j e M 3 : deformation gradient 

V M = | J^- j £ M 3 : displacement gradient 
C = VTV^, e S 3 : right Cauchy-Green strain 



tensor 
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B = V V£ € §2, : left Cauchy-Green strain 
tensor 

E = E{u) = \{C - 1) = ^{Vu T + Vw + Vu T Vu): 
Green-St Venant strain 
tensor 

1 T 

e(u) = 2^ + : nnear i ze d strain tensor 
Pr(Xr) £ R : density per unit mass at 

1>r(Xr) e R 3 : body force density per unit 
mass at Xr € 38 r 
f = PR b R : n -» R 3 
hR(Xg) e R 3 : applied surface force density 
per unit area of 838r at 
X R e d38 R 
g = t XR :Y y ^ R 3 
t R (X R ,n R ),X R € 38 r, \ur\ - 1 : first Piola-Kirchhoff stress 

vector in 38r 
Tr(Xr) : first Piola-Kirchhoff stress 
tensor at Xr € 38 r 
(/, /) - Tr : fi -» M 3 

7# : response functions for 
/ A- = 7 «(/• ). /' e M 3 
^(X*) = V^XrT x T r {X r ) : second Piola-Kirchhoff 

stress tensor at Xr e 38r 
<fij)_= z R : n -» s 3 
Zr> E/?> Z# : response functions for = 
Zr(^) = Z*(0 = TT R (E), 
with F e M 3 , 
C - F r F = I + 2E 
o-*(E) : - 

^(tr£)/ + 2ju£ + 0(£) 
/l,ju : Lame's constants 

v = : Poisson s ratio 

2(A + /u) 
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H(3A + 2fi) 



Young's modulus 



A + fi 

o-ijkt - A5ij5kc + 2fi6jkSje : elasticity coefficients for 

isotropic materials 
W : stored energy function 

-(fi)^ fi (F),FeM 3 + 

^ : stored energy function in 
terms of 

C = F T F(W(F) = ^(O) 
3^* : stored energy function in 
terms of E(W(F) = 
W(I + 2E) = W*(E)) 
(p : stored energy function in 
terms of 

l c : W(F) = cf>(l c ),C = F T F 
: strain energy 

m R n 

Ity) = Wty) - [Btofr) + T\ty)} : total energy 

B(ifj) - J PRbR ■ if/dXR = J f ■ if/dx (for dead loads) 

TlW) ~ J tut' ifrdAR = J g ■ tyda (for dead loads) 



Index 



A 

applied surface force, [5] 

applied surface force. fT8l 

axiom of material frame indifference. 1211 

axiom of static equilibrium.l9l ll9l 

B 

body force, [5] 
body force , fT71 

boundary condition of traction, [54] 
boundary condition of pressure. 1551 
boundary condition of place. l54l 

c 

Cauchy stress vector. [1*31 
Cauchy stress tensor. fT2ll29l 
Cauchy stress vector.l9l l21l 
Cauchy's Theorem. llQI 
Cauchy-Green strain tensors, [3] 
conservative applied surface force. [T%1 
conservative body force. IT8l 
constitutive equation, |20| 

D 

Da Silva's Theorem. 1 191 
dead load.fT71 
deformation.121 1^21 
deformation gradient,|5| 
deformed configuration, [5] 
displacement, |2| 
displacement gradient, [3] 
divergence of a tensor, [5] 

E 

elastic material. l2Tj| 
elasticity coefficient, |36| 



energy. 1^91 

equilibrium equations. I16II37I 
Euler equation. 1391 
Euler's method. I5T1 
Eulerian variable. fT4l 

F 

force: applied surface.l9l l!8l 
body,|9l[T7l 
surface, [5] 
system of,|5| 

G 

Green-St Venant strain tensor, Q 
Green-St Venant strain tensor. l30l 

H 

homogeneous material, |5U| 
homogeneous deformation. 1321 
hyperelastic material. 1351 

I 

incompressible material. l34l 
incompressible material. lll)7l 
incompressible material. lrj2l 
incremental method. [SHI 
infinitesimal rigid displacement. 1191 
isotropic response function. 1261 
isotropic material, [23] 
isotropic response function. 1231 
isotropic stored energy function. 1401 1441 

K 

Korn's inequality. 1641 
L 

Lagrangian variable. fT4l 
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LAME's constants. ETni50irRBl 

linearized strain tensor, [33] 

linearized system of elasticitv. l62M78ll89l 

M 

material 

elastic.l20l 

homogeneous. 1201 

hvperelastic. 1381 

incompressible. l34ll62lfT07l 

isotropic . 1231 

MOONEY-RIVLIN.I9H 

non-homogeneous, 1201 

OGDEN's. l92lfT04l 

St Venant-Kirchhoff, l35ll^ll68llMl 
|g9ll9H[T06l 
material frame indifference. l21ll251l29l 

|40||43l|89l 
mixed displacement-pressure problem, 

E3ED 

mixed displacement-tranction problem, 
l54ll7HllMlfTU3l 

MOONEY-RIVLIN material,^ 

N 

natural state. 1311 
non-homogeneous material. Eol 

o 

OGDEN's material. l92lfT04l 
P 

Piola identitv.l7ll98l 

Piola transform.l4lfT4l 

Piola-Kirchhoff stress tensors. U4IU5ll2^1 

Piola-Kirchhoff stress vector. fT6l 
Poisson's ratio. 1351 
polar factorisation of a matrix. l2Tj| 
polyconvix stored energy function. 1911 

l92lfT00l 
pressure.fT71f!^ll3llfTm 
principal invariants of a matrix. l25l 
principal stretch. |92| 1 1061 
principle of virtual work. fT4| 1 171 1 191 
pure displacement problem. l54M68ll%UI 



pure pressure problem. l55l 
pure traction problem. 1541 

R 

reference configuration, Q 
response function. 1201 
rigid deformation, [6] 
Rivlin-Ericksen Theorem. 1261 

s 

St Venant-Kirchhoff material. l35ll49ll68l 
l84ll89ll93lfT06l 

stored energy function. 1391 
strain 

- energy, [39] 

CAUCHY-GREEN-tensors,[3] 
Green-St Venant-tensor.r7ll30l 
linearized -tensor. l63l 

stress 

CAUCHY - tensor. [HI l29l 
CAUCHY - vector.l9l fL3ll2T1 
PIOLA-KIRCHHOFF - tensors.fT4l 
[Bl[29||43l|45l 

PIOLA-KIRCHHOFF - vector, fT6| 
surface force, [21 
system of forces, [5] 

T 

total energy, [39] 

u 

unilateral boundary condition, [55] 
Y 

YOUNG'S modulus, l35l 



